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ABSTRACT 


A  brief  review  of  the  salient  features  of  the  theoretical  investigation 
of  the  c.  w.  bistatic  inverse  scattering  problem  is  presented.  The  effect  of  changing 
the  origin  of  the  coordinate  system  upon  the  convergent  properties  etc.,  of  the 
spherical  vector  wave  function  representation  of  the  near  scattered  field  and  the 
surface  loci  IS  x  E*  =  0  ,  is  discussed.  It  is  pointed  out  that  a  great  deal 
of  analysis  remains  to  be  done  in  this  area.  The  determination  of  the  surface 
of  the  scattering  body  from  knowledge  of  the  local  total  electric  field  is  given. 
Emphasis  is  placed  upon  the  generalization  of  the  condition  E  x  E*  =  0  as  applied 
to  perfectly-conducting  bodies,  to  scattering  surfaces  characterized  by  the 
impedance  boundary  condition.  Properties  of  the  Matrix  inversion  associated 
with  the  determination  of  the  expansion  coefficients  from  far  field  data  are  dis¬ 
cussed.  Some  numerical  results  are  presented,  and  restrictions  upon  the  choice 
of  aspect  angles  are  deduced. 
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I 

INTRODUCTION 

A  brief  review  of  the  salient  features  of  the  theoretical  investigation 
of  the  c.w.  bistatic  inverse  scattering  problem  is  presented  in  section  II. 

In  section  III  the  effect  of  changing  the  origin  of  the  coordinate  system 
upon  the  convergent  properties  etc.,  of  the  spherical  vector  wave  function 
representation  of  the  near  scattered  field  and  the  surface  loci  E  x  E*  =  0, 
is  considered. 

In  section  IV,  a  discussion  upon  the  determination  of  the  surface  of  the 
scattering  body  from  knowledge  of  the  local  total  electric  field  is  given.  Em¬ 
phasis  is  placed  upon  the  generalization  of  the  condition  E  x  E*  =  0  as  applied 
to  perfectly- conducting  bodies,  to  scattering  surfaces  characterized  by  the 
impedance  boundary  condition. 

In  section  V,  properties  of  the  Matrix  inversion  associated  with  the 
determination  of  the  expansion  coefficients  from  far  field  data,  are  discussed. 
Some  numerical  results  are  presented,  and  restrictions  upon  the  choice  of  aspect 
angles  are  deduced. 
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REVIEW  OF  INVERSE  SCATTERING  THEORY  DEVELOPED  FOR 


C.W.  BISTATIC  SYSTEM 


2.1  THE  PLANE  WAVE  EXPANSION 

It  is  shown  here  that  the  total  field  produced  by  a  plane  wave  Incident 
upon  a  scattering  body,  can  be  expressed  as  the  sum  of  two  terms,  the  incident 
field  and  the  Fourier  transform  of  a  quantity  which  is  related  to  the  scattering 
matrix.  The  resulting  expression  is  valid  for  all  points  in  space  including  the 
interior  of  the  scattering  body. 

To  begin,  the  analysis  will  be  restricted  to  non-magnetic  bodies  (although 
it  could  be  easily  generalized  to  include  such  cases)  and  the  geometry  of  the 
scattering  body  will  be  limited  by  placing  certain  restrictive  analytical  properties 

on  the  surface  S  which  encloses  the  volume  V  of  the  scattering  body.  The 

3 

classes  of  surfaces  chosen,  will  belong  to  class  C  defined  by  Barrar  and 
Dolph  as  follows: 

3 

A  surface  S  is  said  to  belong  to  class  C  if  there  exists  a  finite  number 

m  of  images,  xV  =  x^u,  v),  yV  =  y^u,  v),  z  =  zV (u,  v),  v  -  1.2.3.  .  .  m  of  the 
2  2 

disk  u  +  v  <  1  that  cover  the  surface  S,  and  such  that  the  third  derivatives  of 
x  ,  yU  and  z  with  respect  to  u  and  v  exist  and  are  continuous. 

Harmonic  time  dependence  exp(-iwt)  will  be  taken  in  which  case  Maxwell's 
equations  become 


V  a  E  =  ik  //u  /  e  H 

—  -  o'  o  - 


k^  I  inside  the  body 


where  k^  ,  the  square  of  the  propagation  constant  in  the  body,  is  given  by 


2 


The  incident  electric  intensity  will  be  expressed  in  the  form 


_i  /0  .-3/2  ik  •  x 
E  =  (27 r)  '  a  e  -  - 


(2.1) 


where  k  is  the  direction  of  the  incident  wave,  and  a  ,  the  unit  vector  denoting 
polarization. 

The  homogeneous  body  will  be  treated  first.  In  this  case,  it  follows 
from  Barrar  and  Dolph,  that  the  total  electric  intensity  E  induced  by  the  plane 
wave  Eq.  (2.1),  incident  upon  the  body,  satisfies  the  following  integral  equation 
(which  has  a  unique  solution) 


i  (kl  "  *0 
E  (x)  =  E  (x)  + 


477 


0  E  dv  -  ~  V  •  /  0(E~- n)ds 


V 


(2.2) 


where 


0  = 


ikR 

) 

~R~ 


R  =  (x  -  x') 


and  E  is  the  value  of  E  obtained  by  approaching  the  surface  from  the  interior 
of  the  body.  If  E+  is  the  value  obtained  by  approaching  the  surface  from  the 


exterior  of  the  body,  it  follows  from  the  continuity  relations  that 


2  -  2  + 
k:  E  .  n  =  k  E  •  n 


The  appropriate  expression  for  the  magnetic  field  is  given  by 

... 


(2.3) 


H(x)  =  H  (x)  + 


—  —  477  iu/u 


V  A  I  E  0  dv 


(2.4) 


Before  deriving  a  plane  wave  representation  for  E  (x),  the  following 
vector  will  be  introduced. 
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From  Eq.  (2. 2),  it  is  seen  that  when  |x|  -*  oo  ,  in  the  direction  given  by  the 
vector  k'  ,  then  the  total  field  becomes 

. —  ik  |  x  | 

E(x)  ~  E1(x)  +  /—  ~ — 7-  T(k',k)  (2.6) 

- 1  I  --  v  2  X 

x  I— ►(»  1  ~  1 


indicating  that  the  vector  T  (k',  k)  is  related  to  the  scattering  matrix.  When 
Ik'l  =  |k|  ,  T(k',k)  is  a  measurable  function  being  proportioned  to  the  far 

scattered  field  in  direction  k'  ,  which  is  produced  by  a  plane  wave  of  frequency 

A 

ck  ,  incident  upon  the  body  with  direction  of  incidence  given  by  k  .  The  following 
theorem  may  now  be  proven. 

Theorem: 


If  t  (k1,  k)  = 


1 


(2t r)3/2 


then  <£(k'-k)  =  6(k-k')a-  lim 


-ik'-  v 

e  -  -  E  (x,  k)  d x  , 


T  (k’,k) 


2  2 
e  o  k  -  k'  +  ie 


(2.7) 


(2.8) 


Proof: 


It  follows  from  Eq.  (2.2)  that 


f  (k?  '  k2) 

6  (k1,  k)  =  6(k-k')a+  - — - 


(277) 


3/2 


E  (x',  k)  I  e  —  0dxdv 


-ik'  *x 
e  —  — 


(2tt) 


3/2 


0  (E  •  n)  ds 
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On  setting 


<P 


0  (E  •  n)  ds 


one  can  show  that 


-ik'-  x 

e  -  -  V  ip  dx 


lim 

R  — »  co 
oo 


-ik'-  x  ...... 

e  —  ~  ip  ds  +  ik' 


-ik'-  x 
e  -  ~  p  dx 


oo 


when  S  is  the  surface  of  a  sphere  of  radius  R  .  Letting  k  have  a  small 
oo  oo 

imaginary  part,  it  is  seen  that  the  surface  integral  will  vanish  when  R  oo  . 

oo 

The  resulting  integrals  can  be  reduced  as  follows 


(27r) 


3/2 


-ik'-x  [2  1 


dp  dx 


(2  n)" 


•  (x  -  x')  -  ik'-x 
2  .2 

p  -  k  -  le 


2  ‘d£  6 


-ip  •  X 


n 


2  i2  • 

p  -  k  -  ic 


6  (p-  k') 


.  -ik'  •  x 
2  e  -  - 


*  k2  -  k'2  +  ie 


Combining  the  above  expressions,  one  obtains  the  result 

C  X  (k1  k) 

C_  (k' •  k)  =  <5  ( k  -  k')  a  -  lim  ~  ~  - 

e  -v  o  k  -  k'  +  ie 

As  an  immediate  consequence  of  the  above  theorem,  the  total  field  everywhere 
in  space  including  the  interior  of  the  body, can  be  expressed  in  terms  of  the  incident 
field  and  the  quantity  T(k',  k),  as  follows 


5 


E  (x,  k)  =  E  (x,  k)  + 


(2n) 


3/2 


ip  •  x 

2  i2 
p  -  k 


T  (£,  k)  dj> 


(2.9) 


If  the  integration  space  jd  is  expressed  in  spherical  polar  coordinates  (p,  6  ,  0  ) 

P  P 

where  the  range  of  the  variables  are  -oo<p<oo,O<0  <  7r  /  2  ,  and 

p 

0  <  0^  <  2n  ,  then  in  the  above  integral,  the  contour  of  the  variable  p  bends 
above  the  pole  at  p  =  -k,  and  below  the  pole  at  p  =  k  . 

It  can  be  shown  in  a  similar  manner  that  the  magnetic  field  can  be 
expressed  in  a  similar  form 


H  (x,  k)  =  H*  (x,  k)  + 


u  /u  (2n) 
o 


3/2 


J3  A  T  (jd,  k)  dj3 


(2.10) 


The  above  results  were  derived  for  a  homogeneous  body.  The  results 

may  be  extended  to  include  inhomogeneous  non-conducting  bodies,  i.e.  where 
2  2  2 

kj  =u  e(i  (a  real  function)  and  k^  varies  continuously  in  the  medium.  In  this 
case,  the  appropriate  integral  equation  for  the  total  electric  field  is  (Barrar  and 
Dolph) 


Unfortunately  from  the  standpoint  of  rigor,  only  the  uniqueness  of  the 
solution  of  this  integral  equation  has  been  proven.  It's  existance  has  not  yet 
been  shown. 

Additional  relationships  involving  T  (k1,  k)  have  been  shown  in 
8579-4-Q  ( 1968),  among  which  is  an  integral  equation  involving  T  (k',  k) 
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and  the  Fourier  transform  of  the  body. 

2.2  ALTERNATIVE  REPRESENTATION  FOR  POINTS  EXTERIOR  TO  THE  BODY 
For  a  fixed  direction  of  incidence  and  frequency,  the  quantity  T  (k',  k) 
is  required  for  all  values  of  k'  ,  in  order  to  obtain  expressions  for  the  field 
everywhere  in  space.  Unfortunately,  for  fixed  k,  T  ( k ' , _k)  is  a  measurable 
function  (far  scattered  field),  only  for  values  of  k'  such  that  |  k'  |  =  |  k  |  . 
However  it  will  be  shown  that  restricting  the  requirements  on  the  knowledge 
of  T  (k*»  k)  to  the  values  of  k'  which  lie  on  the  sphere  J  k'  \  =  |  k|,  the  total 
field  can  be  obtained  everywhere  outside  the  minimum  convex  surface  enclosing 
the  body.  To  show  this,  let  the  plane  z  =  zq  be  the  tangent  plane  of  the  body, 
such  that  the  body  lies  in  the  half-space  z  <  z^  .  For  points  in  the  half-space 
z  >  z^,  the  following  representation  of  the  near  scattered  field 


(2.11) 


can  be  reduced  as  follows: 


Set  p  =  k  sin  a  cos  JS 


x 


Py  =  k  sin  a  sin  /3 


Pz 


when  the  domains  of  integration  are 


0  <  ft  <  2n  ,  0  <  a  <  7r /2  -  ioo  ,  and  -  oo  <  q  <  oo 


Expression  (2.11)  becomes 


n  n / 2  -  ioo  r  2n  poo 


E 


s 


k 


q  -  cos  a  -  ie 

Jo  v  -oo 


g (a,  (},  q)  coso'SinQ'do'd^dq, 


(2.12) 
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From  expression  (2.5),  it  can  be  seen  that  for  z  >  zq,  the  contour  in  the 
following  integral 

gikzq  ^  ^  y  dq 
2  2 

,  q  -  cos  a-  ie 

v  -oo 


may  be  deformed,  to  yield  the  following 

.  ik  cos  a  z  _ _  ,  ,  .  „  , 

7r  i  e  T  (k  sin  a  cos  k  sin  a  sin  /3,  k  cos  a ;  k) 

Hence  expression  (2.12)  may  be  placed  in  the  form 

2ir 

e*-  —  Eq  (a,  /3)  sin  a  d  a  d  /3 
o 


1  n  j 2  -  ioo 


E 


ik 

2ir 


where 


(2.13) 


k'  =  k  (sin  a  cos  /3,  sin  a  sin  j3,  cos  a) 

and  Eq  is  the  amplitude  and  phase  of  the  far  scattered  field  in  direction  given 
by  (6,0),  i.e. 

ikR 

ES~  Ar—  E  (0,0).  (2.14) 

irv  O 

Thus  it  can  be  seen  that  for  c.w.  bistatic  scattering,  where  E  (0,0)  is 

— o 

a  measurable  function,  the  scattered  field  for  all  points  in  the  domain  z  >  z  , 

o 

can  be  obtained  from  expression  (2. 13)  provided  that  the  analytic  continuation 
of  E^  (0, 0)  can  be  computed  for  complex  values  of  0  ,  from  knowledge  of  real 
values  of  0  .  It  has  been  shown  (Weston,  Bowman,  Ar,  1966  )  that  the  analytic 
continuation  to  complex  values  of  0  is  possible  provided  that  E  (0,  0)  is 
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known  over  the  complete  domain  0  <  0  <  7r  ,  O<0<27r,  and  a  procedure 
for  doing  so  is  presented. 

By  rotating  the  coordinate  system,  equivalent  expressions  to  (2. 13) 
may  be  derived  which  yield  the  scattered  field  everywhere  outside  a  different 
half-space  excluding  the  body.  Since  the  union  of  all  the  half-spaces  exterior 
to  the  body  is  bounded  by  the  minimum  convex  surface  enclosing  the  body,  it 
can  be  deduced  that  the  scattered  field  at  any  point  outside  this  convex  surface 
can  be  obtained  from  the  knowledge  of  the  far  field  bistatic  data  (measured 
over  all  angles). 

At  this  point,  it  should  be  mentioned  that  expression  (2. 13)  is  not  the 
only  expression  that  can  be  used  to  evaluate  the  near  scattered  field  from  far 
field  data.  It  has  been  shown  that  expression  (2.13)  can  be  re-expressed  in 
terms  of  an  expansion  of  the  spherical  vector  wave  functions.  For  a  fixed 
region  of  the  coordinate  system,  the  expansion  will  be  convergent  down  to  at 
least  the  minimum  sphere  enclosing  the  body.  Alternatively  it  was  shown 
(8579-2-Q)  that  expressions  (2.  13)  could  be  placed  in  the  form 


E 


s 


R  n  En  (0,  0) 


where  E  are  derived  from  the  E  ,  by  a  recurrance  relation  (Wilcox  expansion), 
-n  -n- 1 

2.3  EXTENSION  OF  THE  DOMAIN  OF  CONVERGENCE  AND  THE  CONCEPT  OF 
EQUIVALENT  SOURCES 

It  was  shown  above  that  for  a  fixed  direction  of  incidence  (fixed  transmitter 
position),  one  could  derive  the  near  field  everywhere  outside  the  minimum  convex 
surface  enclosing  the  body,  from  measurements  of  the  phase,  amplitude  and 
polarization  of  the  far  scattered  field  made  over  all  directions.  This  result 
has  been  extended.  Considering  representation  (2. 13)  it  was  shown  that  it  was 
convergent  down  to  the  tangent  plane  z  =  z^,  and  represented  the  scattered 
field  everywhere  on  the  half-space  z  >  zq  .  However  for  perfectly- conducting 
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smooth  convex  bodies,  it  has  been  shown  (Weston,  Bowman, Ar,  1966)  that  this 
expression  was  convergent  down  to  the  plane  z  =  z*  where  z*  >  z^  ,  where 
the  value  of  z*  depends  upon  the  global  geometry  of  the  body.  Since  the  true 
total  field  displays  a  discontinuity  across  the  physical  surface  of  the  body, 
expression  (2. 13)  combined  with  the  incident  field,  cannot  represent  the  true 
field  in  the  intersection  of  the  half-space  z  >  z*  and  the  volume  of  the  body. 
Therefore  in  this  region,  expression  (2. 13)  is  said  to  represent  the  field 
arising  from  a  set  of  equivalent  sources.  In  particular  it  was  shown  that  equivalent 
sources  for  a  prolate  spheroid  lie  on  a  line  joining  the  focal  point,  and  for  a 
sphere,  the  point  at  the  center.  Thus  for  smooth  convex  bodies,  the  domains 
of  convergence  are  not  limited  by  the  minimum  convex  surface  enclosing  the 
body,  but  by  a  much  smaller  volume,  that  is  contained  by  the  minimum  convex 
surface  enclosing  the  equivalent  sources.  As  an  example,  for  the  prolate 
spheroid  with  its  center  lying  at  the  origin  of  the  coordinate  system  and  oriented 
so  that  its  major  axis  lies  along  the  z-axis,  representation  (2. 13)  is  convergent 
in  the  half  space  z  >  c  ,  where  c  is  the  focal  length,  and  both  Wilcox's 
representation  and  the  representation  involving  the  vector  spherical  wave  functions^ 
will  be  convergent  down  to  the  sphere  of  radius  c. 

The  extension  of  the  domain  of  convergence  will  hold  for  piecewise 
smooth  conducting  bodies.  However  the  minimum  convex  surface  enclosing  the 
equivalent  sources  must  enclose  the  curves  or  points  of  surface  discontinuties 
such  as  edges. 

For  cavity  regions  that  penetrate  the  minimum  convex  surface  that  en¬ 
closes  the  equivalent  sources,  alternative  representation  for  the  scattered  or 
near  field  must  be  employed.  Such  representations  will  not  be  of  the  exterior 
type,  the  plane  wave  representation  (2. 13),  the  expansion  in  vector  spherical 
wave  functions,  or  Wilcox's  expansion,  which  are  based  upon  knowledge  of  the 
far  scattered  field  as  a  starting  point,  but  must  be  of  interior  type,  i.e.  based 
upon  knowledge  of  the  total  field  in  a  domain  adjacent  to  the  body.  Such 
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representations  are  only  convergent  in  a  finite  volume.  One  particular  type 
of  such  a  representation  is  given  in  7466-1-F  (Weston  Bowman,  Ar,  1966) . 

2.4  DETERMINATION  OF  THE  SURFACE  OF  PERFECTLY-CONDUCTING 
BODIES. 

Given  expressions  for  the  total  field  (obtained  from  far  field  data) 

which  can  be  computed  in  the  vicinity  of  the  surface  of  the  body,  the  next  step 

is  to  employ  techniques  which  will  locate  the  surface  of  the  body.  For  perfectly 

conducting  bodies,  the  tangential  component  of  E  (the  total  field)  must  vanish. 

This  implies  that  at  a  point  on  the  surface  E  =  n  £  where  n  in  the  unit  normal 

to  the  surface  and  £  is  a  complex  quantity.  For  a  general  point  in  space  the 

total  field  can  be  decomposed  into  the  real  vectors  <£  and  £  as  follows 

I  R 


e  ■ 


For  a  point  on  the  surface  £  a  n  =  0  and  £  a  n  =  0,  which  implies  that 

""I  R 

£-  A  £  t,  =  0  .  Thus  a  necessary  condition  for  a  point  to  be  on  the  surface  of  a 
I  R 

perfectly  conducting  body  is  that,  at  that  point,  the  total  field  must  satisfy  the 
condition 


E  A  E*  =  0  , 


when  E*  is  the  complex  conjugate  of  E  .  The  advantage  of  this  condition  is 
that  it  is  a  local  condition,  requiring  only  the  calculation  of  the  total  field  at 
the  point  in  question.  Since  the  above  condition  is  not  sufficient,  additional 
requirements  would  be  required,  such  that  E  x  E*  must  vanish  at  adjacent  points 
x  and  x  +  Ax  ,  and  that 


lim 

A  x  0 


E  (x)  •  Ax 

“FIT 


0 


which  implies  that  E  is  normal  to  surface.  In  addition,  the  resulting  surface 
formed  by  the  set  of  points  must  be  closed.  However,  even  if  a  surface  S 

o 
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was  found  such  that  it  satisfied  these  conditions,  it  would  not  necessarily  be 
the  correct  surface.  Since  for  the  enclosed  volume  formed  by  two  closed 
surface,  on  which  n  x  E  vanishes,  there  is  a  discrete  spectrum  (in  frequency), 
this  implies  that  at  a  particular  frequency,  there  may  exist  additional  surfaces 
for  which  n  x  E  vanishes.  However  the  geometry  of  the  additional  surfaces 
depend  upon  frequency,  and  these  can  be  separated  out  from  the  proper  surface 
by  employing  at  least  two  different  frequencies. 

For  the  illuminated  region  of  the  body,  an  auxiliary  condition  was  de¬ 
veloped  which  considerably  helped  discriminate  the  proper  surface  from  the 
loci  E  x  E*  =  0.  This  condition  given  as  J  E*  j  -  |  ES  |  =  0  ,  yields  an  approximate 
surface,  which  approaches  smooth  convex  portions  of  the  correct  surface  in 
the  limit  of  high  frequency  scattering. 

From  the  numerical  standpoint  where  there  are  errors  due  to  input 
information,  etc.,  it  was  demonstrated  both  theoretically  and  numerically 
that  the  condition  E  x  E *  =  0,  should  be  replaced  by  finding  the  minima  of 
|E  x  E*)  . 

2.5  FAR  FIELD  INFORMATION  LIMITED  TO  A  SOLID  ANGLE. 

For  most  practical  solutions,  the  far  scattered  field  (phase,  amplitude 
and  polarization)  is  measurable  only  over  a  region  of  limited  bistatic  angles. 

In  this  case  it  is  important  to  know  what  information  can  be  determined  about 
the  body,  where  measurements  of  the  scattered  far  field  are  confined  to  a 
limited  domain  of  aspect  angles.  This  in  turn  depends  upon  the  accuracy  to 
which  the  near  scattered  field  can  be  computed  using  the  limited  far  field  data. 

A  general  discussion  of  this  latter  point  is  presented  in  8579- 1-Q.  However, 
the  employment  of  high  frequency  asymptotic  results  yields  a  far  more  illuminating 
picture  (8579-2-Q).  These  results  will  be  briefly  presented  here.  When  the 
scattered  far  field  is  measured  over  a  limited  aspect  region  (example:  the  cone 
0  <  0  <  0q  ),  the  analytic  continuation  of  quantity  Eq  (a,  fi)  given  in  expression 
(2.13),  cannot  be  found  for  the  complete  complex  a  plane.  In  this  case  the 
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following  approximation  to  expression  (2. 13)  will  be  employed 


E  =  — 


ik 
2t r 


1  277 


0 


6 


J  0 


i  kf  •  x 

e  —  —  E  (o,  6)  sin  a  d  a  d  ft  .  (2.15) 

—  o 


For  the  case  of  the  perfectly  conducting  sphere  of  radius  a  ,  illuminated  by 
the  plane  wave 

„i  a  -ikz 
E  =  x  e 

the  geometric  optics  scattered  far  field  has  the  form 

E  ^  (a,  /3 )  =  -  e  (a,/3)  exp  £-2ika  cos  (0/2)]  , 

where 


e  (o,  j3)  =  x  |[cos  a  cos^j3  +  sin^J  -  y  jj.  -  cos  o]  sin  j3  cos  /3  - 
-  z  sin  a  cos  0  . 

From  relation  (2. 15)  the  near  scattered  field  is  given  approximately  by 

f  2tt  re 


ES(x)  = 

—  ~  47T 


ikf  (a,  ft)  a  .  .  .  ,  , 

e  e(a,p)  smo'do'dfi  , 


(2.16) 


0  J  0 


where 


f(a,j3)  =  r  £sin  6  sin  a  cos  (0-  8)  +  cos  6  cos  oj  -  2a  cos  ^  . 

As  k-^00,  the  dominant  contribution  to  the  integral  arises  from  the  vicinity 
of  the  stationary  phase  point  (j3=0,  a  =  a  )  where  satisfies  the  equation 


r  sin  ( a  -6)  =  a  sin  (tt) 

O  x  Z  s 
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provided  that  0  <  a  <  6  .  By  means  of  first  order  stationary  phase  evaluation 

—  o  “  o 

we  obtain  immediately 


fi 

ES(x)  = 


D(o) 

D(s) 


1/2 


iks  -  ika  cos 


e  (a  ,  0 )  e 
—  o 


(f) 


where 


D(s)  =  (cos  f  +  7-)(l+  fcos  y) 


The  distance  s  given  by  the  relation 


s 


r  cos  (a 

o 


-  0)  -  a  cos 


is  the  distance  along  the  reflected  ray  from  the  point  of  reflection.  The  resulting 

A  g 

expression  for  E  (x)  is  the  geometric  optics  near  field  expression.  When  the 
far  field  information  is  limited  to  the  cone  0  <  0  <  0^  ,  the  near  field  expression 
will  be  limited  to  a  volume  of  space  such  that 


r  sin  (a  -  0)  =  a  sin  ( TT  )  0  <  a  <  0  ,  (2.17) 

where  (r,  0, 0)  are  the  coordinates  of  a  point  in  this  volume.  In  order  to  find 
the  portion  of  the  scattering  surface  which  can  be  determined,  set  r  =  a,  in  which 
case  it  is  seen  that 


0  = 


0  <  a  <  0 
—  o  —  o 


This  then  implies  that  the  near  scattered  field  given  by  expression  (2. 15)  can  be 
found  on  the  portion  of  the  sphere  given  by  the  cone  0  <  0  <  0^  /2  .  For  points 
outside  the  volume  of  space  given  by  Eq.  (2. 17),  stationary  phase  techniques 
cannot  be  employed,  in  which  case  the  asymptotic  approximation  to  Eq.  (2. 16)  is 
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obtained  by  alternative  means . 

The  above  example  illustrates  the  fact  that  in  the  case  of  high  frequency, 
asymptotic  evaluation  of  expression  (2. 15)  is  equivalent  to  tracing  back  the  rays 
to  the  portion  of  the  scattering  surface  from  which  they  arose.  The  result 
is  not  confined  to  just  perfectly  smooth  convex  shapes  but  can  be  applied  to  piece- 
wise  smooth  convex  shapes,  the  far  scattered  field  is  decomposable  into  components 
of  the  form 

Eik  ip  (0,  0) 

A  n(0,0)  e 
n 

and  stationary  phase  evaluation  of  the  expression  (2. 15)  for  each  individual 
component,  effectively  traces  back  the  various  rays  to  their  originating  portions 
of  the  scattering  object.  The  result  is  applicable  to  flat  portions,  in  which 
case  the  far  field  component  approaches  a  delta  function  of  the  angular  variables 
as  k  ->  oo  .  The  example  of  the  flat  plate  given  in  8579-2-Q  illustrates  this. 
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Ill 


EFFECTS  OF  A  DISPLACED  ORIGIN  ON  CONVERGENCE  PROPERTIES 


In  employing  the  Wilcox  theorem  or  the  vector  spherical  wave  function  expansion, 
it  can  be  shown  that  the  near  field  representation,  derived  from  the  far  scattered 
radiation  pattern,  is  uniformly  convergent  for  all  points  outside  the  minimum 
sphere  enclosing  the  equivalent  sources  of  a  smooth  convex  shaped  body.  From 
the  knowledge  of  the  near  field  everywhere  outside  this  minimum  sphere,  it  is 
sought  to  find  the  associated  scattering  surface  by  proper  application  of  suitable 
boundary  conditions.  These  Inverse  Scattering  Boundary  conditions,  in  general, 
ought  not  involve  any  constituent  parameters  (surface  normal  n  ,  surface  locus 
S(x,  y,  z),  relative  impedance  rj  )  of  the  unknown  scatterer  and  must  be  given 
solely  in  terms  of  the  known  nearfield  representations  of  the  incident  and  the 
scattered  field  with  respect  to  a  fixed  origin  which  lies  within  the  scatterer. 

For  a  perfectly  conducting  surface  it  was  found  that  the  following  two  boundary 
conditions  can  be  successfully  applied: 

{l  E1!-  |ES|)  =  0  (3.1) 

i  s 

where  E  denotes  the  incident  field  and  E  the  nearfield  representation  of 
the  scattered  field.  This  boundary  condition  represents  the  geometrical  optics 
limiting  approximation  and  in  general,  is  applicable  only  within  a  narrow  cone 
about  the  specular  point  (8579-3-Q,  8579-4-Q). 

k)  (  t  T* ")  T  i  s 

(E  x  E  j=0,  where  E  =  E+E  (3.2) 

This  boundary  condition  was  derived  by  V.  H.  Weston  (Final  Report  7644-  1-F, 
Weston,  Bowman,  Ar).  Although  it  is  only  a  necessary  but  not  sufficient  condi¬ 
tion,  namely  producing  a  family  of  concentric  surface  loci  in  addition  to  the 
proper  one,  its  application  together  with  1 1 EX|  -  |e  sl)  =  0  proved  to  be 
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indispensable  since  that  portion  of  the  proper  locus  within  the  conical  section 
exterior  to  the  Wilcox  minimum  sphere  can  be  determined  with  great  accuracy 
for  a  minimum  number  of  expansion  terms,  (8579-3-Q,  8579-4-Q). 

To  be  more  precise,  the  combined  boundary  conditions  will  yield  the 
proper  surface  locus  only  on  those  portions  of  the  surface  for  which  the  chosen 
expansion  for  the  scattered  field  is  convergent,  for  a  prolate  spheroid,  the 
boundary  conditions  applied  to  the  complete  vector  spherical  wave  function  ex¬ 
pansion  associated  with  the  origin  located  at  the  center  of  the  body,  will  yield 
portions  of  the  surface  outside  the  sphere  r  =  c  where  c  is  the  semi-focal 
length.  In  order  to  obtain  the  side  portions  one  has  to  displace  the  origin  in 
the  plane  passing  through  the  center  of  the  body  and  perpendicular  to  the  axis  of 
the  body.  Thus  in  using  far  field  data  one  needs  to  know  the  domain  of  convergence 
of  the  expansion.  Generally  this  can  be  only  obtained  if  the  complete  expansion 
were  known,  which  in  practice  is  not  achievable  since  the  number  of  receiver 
locations  will  be  a  finite  number.  If  the  domain  of  convergence  is  not  known,  then 
the  application  of  the  boundary  conditions  to  the  expansion  may  yield  some  correct 
portions  of  the  body,  but  what  the  correct  portions  are  cannot  be  prescribed  unless 
some  additional  condition  is  employed.  However  if  it  was  known  apriori  that  the 
origin  of  the  coordinate  system  was  located  at  the  center  of  the  body,  and  the  portions 
of  the  body  farthest  removed  from  the  origin  where  smooth  and  convex,  then  the 
application  of  the  boundary  condition  will  yield  the  correct  surface  locus  (to  within 
some  numberical  error)  for  these  portions.  Thus  if  the  approximate  location  of 
the  origin  was  not  known,  some  additional  criteria  should  be  employed,  such  as 
the  one  described  below. 

A  possible  additional  condition  would  be  to  displace  the  origin  of  the 
coordinate  system,  and  on  employing  the  boundary  condition  to  the  scattered  field 
expansion  derived  in  the  displaced  coordinate  system  to  obtain  a  new  surface 
locus,  one  would  then  determine  the  shift  in  the  surface  locus  due  to  the  displace¬ 
ment.  If  the  shift  is  negligible  for  portions  of  the  surface,  then  those  portions 
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may  be  the  proper  surface.  For  example,  if  the  origin  of  the  coordinate 
system  located  at  the  center  of  a  prolate  spheroid,  was  shifted  slightly 
down  the  axis,  the  loci  for  the  top  and  bottom  of  the  prolate  spheroid  derived 
from  the  boundary  conditions  applied  to  the  near  field  expansions,  should  be 
the  same  as  before  the  shift.  These  remarks  are  strictly  qualitative.  A 
quantitative  study  should  be  made  and  is  outlined  below. 

Let  x  represent  the  points  in  a  cartesian  frame  of  reference  the  origin 
of  which  is  located  in  the  vicinity  of  the  scattering  body  ( its  precise  location 
was  not  established)  Let  a  plane  wave  be  incident  upon  the  body  producing  a 
scattered  field,  the  far  field 

ikR 

ES  -  ±ir  Eo  (0,  0) 

being  measured  at  a  finite  number  (N)  locations  0nj  •  Let  ?S(x,N) 

be  a  finite  expansion  of  the  near  field  derived  from  the  far  field  data  (the  ex¬ 
pansion  may  in  vector  spherical  harmonics,  or  a  Wilcox  type  expansion) 

Suppose  that  the  origin  is  displaced  a  distance  d,  such  that  if  x' 
represents  a  point  in  the  displaced  coordinate  system  then 

x'  =  x  +  d 

In  the  displaced  coordinate  system,  the  far  field  directions  ^0^,  0n|  remain 
the  same,  (i.e.  for  R' — ►oo,  6'  =  9,  0'  -  0),  however  the  far  field  quantity 
E'  (0',  0')  is  related  to  E  (0,  0)  by  the  relation 

E'  (0',  0')  =  E  (0,  0)  exp  (  -  ik  d) 

_  o  —  ■*-  — o  —  — 

fyj  g 

where  k  is  a  vector  in  the  radial  direction.  Let  E'  (x',  N)  be  a  finite  ex¬ 
pansion  of  the  near  field  derived  from  far  field  data  (N  receivers)  in  the  displaced 
coordinate  system. 
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If  x  is  a  point  in  space  where  the  complete  (N  =  oo)  expansions  are 
convergent  the  difference 

ES  (x,  oo)  -  E'  (x  +  d,  oo) 

would  be  zero.  However  when  a  finite  number  of  terms  are  used  in  the  expansions 
(finite  number  of  receivers)  the  difference  will  not  be  zero,  except  of  course  at 
the  receiver  positions  in  the  far  field.  Thus  it  would  be  fruitful  to  determine 
quantitatively  a  number  6  (d)  such  that  if 


lim 

d  — ►  o 


ES  (x)  -  E,S  (x  +  d) 


<  6  (d)  I  E1 


where  E1  is  the  incident  field,  the  application  of  the  boundary  conditions 
on  the  total  field  at  points  x  for  which  the  above  inequality  holds  will  yield  the 
correct  surface  to  a  specified  degree  of  accuracy. 

Such  a  condition  given  above  may  have  to  be  obtained  through  numerical 
examples.  It  will  depend  upon  what  expansion  for  the  near  field  is  chosen,  and 
how  the  expansion  is  matched  to  far  field  data  (i.  e. ,  by  matrix  inversion,  least 
squares,  or  otherwise). 

To  obtain  the  portions  of  the  surface  outside  the  domain  of  convergence 
(if  the  complete  expansion  were  used  or  outside  the  domain  given  by  the  above 
inequality  for  a  finite  expansion,  the  origin  of  the  coordinate  system  would  have 
to  be  significantly  changed.  For  this  new  location  of  the  origin,  the  domain 
validity  of  a  finite  expansion  would  have  to  be  chosen  by  displacing  slightly  the 
origin  and  employing  the  above  inequality. 

Some  initial  investigations  of  the  effects  of  displacing  the  origin  for  a 
finite  expansion  in  vector  spherical  wave  functions,  upon  the  surfaces  derived 

jJc 

through  the  boundary  condition  ExE  =0,  is  discussed  in  internal  memo 
(W.  Boerner).  However,  this  memo  is  just  an  initial  attempt  to  understand  the 
problem  and  will  not  be  reported  here. 
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IV 


ASPECTS  ON  CONTINUOUS  WAVE  INVERSE  SCATTERING  BOUNDARY 
CONDITIONS  FOR  THE  DETERMINATION  OF  THE 
SHAPE  AND  THE  MATERIAL  CONSTITUENTS  OF  A  CONDUCTING  BODY 


4.1  INTRODUCTION 

Although  methods  have  been  outlined  on  how  to  determine  the  material 
characteristics  of  a  conducting  body,  employing  the  monostatic-bistatic  theorem 
(7644-1-F,  Chapter  XI),  no  suitable  method  has  been  found  so  far  which  will 
simultaneously  determine  the  shape  and  the  associated  material  constituents 
of  an  unknown  body,  since  the  inverse  scattering  boundary  condition 
1 1  E  xE  =0,  where  E  =E  +E  j  ,  derived  by  Weston  (7644-1-F, 


Chapter  X)  holds  only  for  the  case  of  a  perfectly  conducting  body.  This  condi- 


though  necessary  but  not  sufficient,  may  however  be 


tion 


considered  as  a  first  step  of  formulating  a  more  general  c.w.  Inverse  Scattering 
Boundary  Condition. 

An  attempt  of  generalizing  this  condition  will  be  prescribed  in  this  section. 
Before  any  C.  W.  I.  S.  boundary  condition  can  be  applied,  a  suitable  near  field 
representation  of  the  scattered  field  from  the  measured  far  field  radar  data 
must  be  sought.  This  can  be  achieved  by  either  employing  a  series  expansion 
into  proper  vector  wave  functions  and  if  so  required  their  associated  plane  wave 
integral  representation  (7644-1-F,  Chapter  V;  8579-1-Q,  8579-2-Q)  or  an 
expansion  method  derived  by  Wilcox  ( 1956)  and  Muller  ( 1956)  as  discussed  in 
7466-1-F  (Chapter  IV).  The  associated  expansion  coefficients  may  be  obtained 
from  a  matrix  inversion  technique  (8579-2-Q  and  8579-1-F).  Assuming  that 
such  a  sufficiently  accurate  nearfield  representation  of  both  the  electric  and  the 
magnetic  field  vector  is  found,  the  question  of  how  to  derive  suitable  I.S.  boundary 
conditions  arises  which  as  well  can  be  applied  to  the  determination  of  the  character¬ 
istic  parameters  of  conducting  bodies  and  their  shapes . 

In  contrast  to  problems  of  direct  scattering  and  diffraction  for  which 
the  shape  and  the  material  constituents  of  the  scatterer  are  assumed  to  be  known 
a  priori  together  with  the  prespecified  incident  field  vector  E*  as  regards  the 
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computational  coordinate  system  and  thus  may  be  incorporated  into  the  boundary 
conditions,  in  problems  of  inverse  scattering  such  boundary  conditions  must  be 
sought  which  in  particular  do  not  depend  upon  either  the  shape  or  the  material 
properties  of  the  scattering  body  and  its  enclosing  surface,  but  allow  to  specify 
those  characteristic  parameters  uniquely.  And  this  solely  from  the  near  field 
representation  of  the  electric  and  the  magnetic  field  vectors,  derived  from  the 
measured  far  field  data  as  described  above.  Now  the  question  arises  as  to 
how  many  and  which  characteristic  parameters  must  be  defined  to  uniquely  deter¬ 
mine  the  shape  and  the  material  constituents  of  the  unknown  scatterer.  In  the 
scalar  case  it  is  sufficient  to  employ  the  following  parameters: 

i)  the  local  surface  normal  n  (R,  0,  0)  of  the  proper  surface, 
ii)  a  relative  surface  impedance  q  (R,  0, 0)  which  is  a  scalar  quantity, 
or  the  interior  propagation  constant  k.^  (R,  0,  0). 

Thus  at  least  three  independent  characteristic  equations,  solely  expressed  in 

T  i  s  T  1  s 

terms  of  the  near  field  quantities  E  =  E  +  E  and  H  =  H  +  H  ,  must  be 
found  to  determine  the  surface  locus  S  (R,  0,  0)  and  the  surface  impedance 
n  =  Ini  exp  i  ^  . 

Whereas  in  the  vector  case  the  relative  surface  impedance  must  be 
formulated  in  dyadic  or  tensor  formulation,  involving  a  further  unknown  parameter 

A  A 

which  may  be  described  by  a  local  polarization  angle  e  (R,  0 , 0;  k.  ,  e  )  or 

A  a  1  * 

a  local  depolarization  angle  6  (R,  0,  0;  k^ ,  e^).  These  angles  can  be  expressed 

in  terms  of  the  surface  normal  n  (R,  0J  0)  as  well  as  the  unit  direction  vector 
•  a  A 

of  the  incident  wave  k^  and  its  associated  polarization  vector  e^  .  Thus  the 

impedance  dyadic  17  will  be  a  function  of  both  the  material  properties  of  the 

scattering  surface,  the  surface  locus  S(Ri,0>0)  or  its  associated  local  normal 

A  A  A 

n  (R,  0,  0),  and  the  prespecified  properties  k^,  e^  of  the  incident  wave.  This 
additional  complication  may  ask  for  one  or  two  more  independent  characteristic 
equations  for  the  elimination  of  e  or  i. 
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If  such  a  set  of  independent  scalar  and  vector  equations  exist  which  can 

be  employed  to  uniquely  determine  the  surface  locus  S  (R,  0,  0)  and  the  material 

characteristic  rj  ,  then  one  may  argue  that  the  I.S.  boundary  condition 
X  x* 

E  x  E  =0  may  constitute  the  remaining  part  of  such  a  set  of  independent 

equations  for  the  degenerate  case  of  r/  =  0  .  With  this  goal  in  mind  more  de- 

r  X  T*  T  i  s') 

tailed  properties  of  <  E  xE  =  0,  E  xE+Er  will  be  discussed  next, 

as  well  as  properties  of  the  Leontovich  boundary  condition  and  its  complementary 

formulations.  Then  it  will  be  shown  that  for  the  scalar  case  a  set  of  independent 

scalar  and  vector  equations  can  be  derived,  based  upon  the  concept  of  describing 

the  material  surface  properties  by  a  Leontovich  boundary  condition,  which  will 

X  X* 

in  the  limit  as  rj  =  0  ,  degenerate  into  the  condition  E  x  E  =  0  as  derived 
by  Weston  in  7644-1-F  (Chapter  X). 

rp  npjV 

4.2  THE  C.W.  I.S.  BOUNDARY  CONDITION  E  xE  =0  FOR  THE  DEGENERATE 
CASE  OF  A  PERFECTLY  CONDUCTING  SCATTERER,  i.e.  r)  =  0  . 

In  chapter  X  of  7644-1-F,  the  necessary  but  not  sufficient  condition 
X  X* 

E  x  E  =0  was  derived  intuitively  based  upon  the  physical  properties  that  the 

T  i  s  a  T 

total  field  E  =  E  +  E  must  satisfy  the  boundary  condition  n  x  E  =  0  on 

A 

the  surface  of  a  perfectly  conducting  body,  and  since  the  surface  normal  n  is 

T  T  T* 

a  real  vector  also  E  cannot  be  of  complex  direction,  thus  E  x  E  =0. 

The  application  of  this  boundary  condition  is  verified  for  spherical  and  prolate 

spheroidal  test  cases  in  8579-3-Q  and  8579-4-Q,  together  with  the  condition 

{IS*  I  -  I  E3/  ]  =  0  as  derived  from  the  geometrical  optics  approximation.  Both 

conditions  yield  the  most  accurate  results  if  the  incident  polarization  is  parallel 

i  s 

to  the  generators  of  the  scattering  body,  since  then  E  +  E  =0  identically. 

X  X* 

Although  E  x  E  =  0  is  only  a  necessary  condition  but  not  sufficient  in  that 
an  infinite  set  of  exterior  concentrical  (hyperbolic)  surface  loci,  and  for  scatterers 
of  larger  electrical  measure  ka  a  limited  set  of  interior  pseudo  loci  (emanating 
from  the  associated  interior  caustic)  in  addition  to  the  proper  one  are  generated, 
this  proper  surface  locus  can  be  determined  with  least  error.  In  fact,  in  8579-3-Q 
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T  x* 

and  8579-4-Q,  it  is  verified  that  with  E  x  E  =0  the  proper  surface  locus 

can  be  identified  more  accurately  as  with  the  condition  1 E 1 1  -  |  E S  |  =  0  for 

which  only  one  locus  is  obtained,  however  with  pronounced  deviations  from  the 

X  x* 

proper  one  (8579-4-Q).  The  family  of  resulting  surface  loci  of  E  xE  =0 

are  equidistantly  spaced,  concentrical  hyperboloids  with  AX  =  kAR  =  n / 2 

spacing  along  the  axis  of  symmetry  for  rotationally  symmetric  bodies  with 

end-on  incidence.  The  resulting  plots  for  one  and  the  same  scatterer  versus 

the  electric  length  are  identical  for  different  frequencies.  If  however,  the 

resulting  loci  are  plotted  versus  the  geometrical  length  R  =  X/k  ,  then  it  can 

be  shown  (8579-4-Q)  that  the  proper  surface  locus  is  stationary  whereas  all  the 

other  interior  as  well  as  exterior  additional  loci  will  shift,  since  the  proper 

surface  locus  is  independent  of  the  wavelength  and  all  further  loci  correspond 

to  a  discrete  set  of  eigenfrequencies  (7644-1-F,  chapter  X).  Thus  the  application 

X  x* 

of  the  boundary  condition  E  x  E  -  0  at  only  two  different  frequencies  may  be 

sufficient  to  uniquely  determine  the  shape  of  a  smooth  convex- shaped,  perfectly 

conducting  body  in  those  regions  which  are  exterior  to  the  Wilcox  minimum 

sphere.  Those  portions  of  the  scattering  surface  within  the  Wilcox  minimum 

X  x* 

sphere  can  be  determined  with  E  x  E  =  0  as  well,  but  not  with  the  condition 
[IS*/  ‘IfilJ-O  ,  if  properties  of  displacing  the  computational  origin  from  the 
geometrical  center  are  employed  at  two  different  frequencies  by  matching  the 
geometrical  length  of  the  displacement  vector  so  that  a  stationary  surface  locus 
within  the  angular  domain  along  the  direction  of  the  displacement  vector  is 
found,  (8579- 1-F  ). 

*p  'p* 

Thus  the  boundary  condition  E  x  E  =0  must  be  considered  indispensable 
in  computational  problems  of  C.  W.  Inverse  Scattering  as  applied  to  the  identi¬ 
fication  of  smooth,  convex- shaped,  perfectly  conducting  bodies. 

T  t* 

Finally  the  question  arises  whether  the  boundary  condition  E  x  E  =0 
may  not  bear  more  physical  information  and  whether  its  generalization  to  conducting 
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bodies  may  reveal  a  relationship  with  Maxwell's  Stress-Energy  Tensor  at  the 
bounding  surface  of  a  scatterer.  It  was  therefore  attempted  to  employ  the  Leontovich 
boundary  condition  of  the  scalar  case  to  study  this  question  in  more  detail. 


4  .3  THE  LEONTOVICH  BOUNDARY  CONDITION  AND  ITS  VARIOUS 
COMPLEMENTARY  FORMULATIONS. 

In  problems  of  C.  W.  inverse  scattering,  the  concept  of  an  impedance 
boundary  condition  may  be  employed  to  its  best,  if  such  an  impedance  can  be 
defined  so  that  it  describes  the  averaged  local  electromagnetic  properties  of 
the  unknown  scattering  surface.  The  boundary  condition  (4.3. 1),  known  as 
Leontovich  condition  (Leontovich,  1948)  suggests  to  offer  the  desired  formula¬ 
tion  for  the  scalar  case: 


where 


ET  -  (n  •  ET)  n  =  n  Zq  n  x  HT 


(4.3.1) 


et  = 


E*  +  ES 


the  total  electric  field 


T  i 
H  =  H  + 


H 


the  total  magnetic  field  in  the  region  surrounding 
the  body. 


— —  =  120  ir  Q. ,  the  intrinsic  impedance  of 
o  free  space. 


n  =  the  unit  outward  normal  as  regards  the  surface  which 
is  a  purely  real  vector. 


r)  ,  the  relative  impedance,  a  complex  scalar  quantity,  designated 
as  Leontovich  impedance. 


The  degenerate  cases  for  which  (4  3.1)  is  not  applicable  are  treated  in  detail 
in  Senior  (1959,  1962)  and  Weston  (CAA-0020-10-TR)  and  will  be  excluded  from 
further  discussion.  The  direct  application  of  (4.  3. 1)  would  imply  the  a  priori 
knowledge  of  n  =  n  (R,  0,  0)  and  r)  =  r)  (R,  Q,  0).  Thus  other  formulations  of 


24 


(4.3.1)  are  required  which,  however,  must  contain  identical  information  such 

that  the  surface  locus  S  (R,  9, 0)  and  q  (R,  9, 0)  can  be  found  solely  from 
T  T 

E  and  H  .  The  first  additional  equation  is  found  by  applying  a  vector  product 
operation  of  n  onto  (4.3.1),  yielding 

E^xn  =  rj  Zq  |jlT  -  (n  •  H^)  nj  (4.3.2) 

which  merely  represents  an  alternative  statement  of  (4.3. 1)  and  may  be  rewritten 
with  5  =  l/rj  into 

HT  -  (n  •  HT)  n  =  -  ?Y  (n  x  ET)  (4.3.2a) 

—  —  O  — 

T  T 

and  thus  corresponds  to  (4.3.1)  under  the  transformation  E  -»  H  , 

T  T 

Z  H  -¥  -  Y  E  )  and  n  5  ,  where  5  denotes  the  relative  admittance, 
o  —  o  — 

Senior  (1962a)  has  shown  the  affinities  of  this  transformation  with  Babinet's 
principle  and  has  proven  the  invariance  of  this  transformation,  attributing  rj 
to  the  material  constituents  of  the  body. 

In  addition  to  (4.3.1)  and  (4.3.2)  its  conjugated  formulations  will  be 
introduced  as: 

T*  a  T*  a  a  T* 

E  -  (n  •  E  )  n  =  rf*  Z  n  x  H  (4.3.3) 

—  —  o  — 

and 

ET*x  n  =  rf  Zq  ^  HT>,t  -  (n  •  HT*)  n]  (4.3.4) 

or 

[^-(n.H^  =  -C*Y  (nxE1^)  .  (4.3.4a) 

L—  —  O  — 

The  validity  of  statements  (4.3.3)  and  (4.3.4)  must  strictly  assume 

T  T 

that  all  implied  field  quantities  E  and  H  ,  as  well  as  the  relative  impedance 
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r)  =  rj  (R,  0,  0)  are  analytical  functions,  and  n  =  n  (R,  0,  0)  is  piecewise  continuous, 

satisfying  the  set  of  linear  equations  (4.3. 1)  and  (4.3.2)  which  in  turn  satisfy 

Maxwell's  equations.  The  particular  constraints  on  q  and  n  for  which  this  is 

not  the  case  are  discussed  by  Weston  (CAA-0020-10-TR),  and  scatterers  of  these 

particular  features  will  be  excluded.  The  conjugation  of  (4 . 3 . 1)  and  (4. 3 . 2) 

solely  implies  the  reversal  of  the  reactive  character  of  all  implied  electromagnetic 

quantities,  thus  the  pair  (4.3.3)  -  (4.3.4a)  bears  similar  affinities  as  does  the 

pair  (4.3.1)  -  (4.3.2a),  now  however,  for  a  surface  of  reversed  reactive  character. 

Assuming  that  the  electromagnetic  behavior  in  the  vicinity  of  a  scatterer 

satisfies  the  Leontovich  boundary  condition  (4.3. 1),  the  question  must  be 

answered  whether  a  set  of  scalar  and  vector  equations  can  be  derived  from  the 

complementary  set  (4.3. 1)  to  (4.3.4)  of  the  Leontovich  boundary  condition 

which  may  be  employed  to  determine  the  surface  locus  S  (R,  0, 0),  the  surface 

normal  &  (R,  0,  0),  the  relative  surface  impedance  rj  (R,  0,  0)  solely  from  the 

T  T 

knowledge  of  the  near  field  representations  of  E  and  H  .  This  may  be 
achieved  by  proper  scalar  and  vector  product  operations  of  (4.3. 1)  to  (4.3.4) 
onto  one  another  which  must  result  in  the  following  set  of  independent  equations: 

i)  A  vector  or  scalar  equation  which  ought  to  be  independent  of  ri  and 
may  thus  determine  the  surface  locus  S  (R,  0,  0)  of  the  unknown  target. 
This  equation  may  possess  a  discrete  set  of  additional  solutions, 
resulting  in  a  family  of  concentric  surface  loci, 

ii)  A  set  of  independent  scalar  equations  necessary  to  define  both  the 
amplitude  /  ri  |  and  the  phase  q//n/=  exp  i  ip  of  the  relative  surface 
impedance.  These  equations  should  result  in  a  pair  of  quadratic  equations 
such  that  a  fourfold  complementary  solution  is  obtained,  corresponding 
to  the  complementary  character  of  equations  (4.3.  1)  to  (4.3.4). 

iii)  A  closed  form  expression  for  the  surface  normal  n  ,  which  must 
yield  a  purely  real  vector  quantity  and  ought  not  depend  upon  ri  , 
otherwise  the  correct  solution  of  rj ,  satisfying  (4.3.1)  cannot  be  found. 
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To  this  end  the  scalar  formulation  of  the  Leontovich  boundary  condition 
was  considered  only  which  in  general  cannot  be  employed  for  problems  of  inverse 
scattering  as  applied  to  the  determination  of  the  shape  and  material  constituents 
of  closed  convex-shaped  bodies.  The  vector  wave  nature  of  the  inverse  scattering 
problem  requires  the  definition  of  a  dyadic  formulation  of  the  relative  impedance, 
where  the  following  equations  must  be  satisfied  on  the  surface  of  a  non-perfectly 
conducting  body 


ET  -  (n  •  ET)  n  =  Z  *n  •  (5  x  HT) 

—  —  o  — 

T  a  ,  T  .a  T.  a  . 

E  xn  =  Z  n  •  (H  -  (n  •  H  )  n  ) 

—  o  — 


(4.3.5a) 


(4.3.5b) 


The  impedance  dyadic  is  symmetric  and  will  contain  all  terms  except  for  the 
case  of  plane  wave  incidence  onto  a  planar  surface  with  invariance  along  one  axis 
normal  to  the  plane  of  incidence  for  which  the  diagonal  terms  are  involved  only. 

In  all  other  cases  a  local  angle  of  polarization  as  well  as  a  local  angle  of  de¬ 
polarization  together  with  the  radii  of  curvature  will  be  involved.  These  prop¬ 
erties  however  may  not  allow  the  formulation  of  the  required  independent  equations 
since  the  impedance  dyadic  will  become  a  function  of  the  radiant  dependence 
along  the  outward  normal  (Morse  and  Feshbach,  1953).  Before  any  further 
assumptions  can  be  made,  the  scalar  case  will  be  investigated  in  detail. 

4.4  FORMULATION  OF  A  SET  OF  INVERSE  SCATTERING  BOUNDARY  CONDITIONS 
FOR  A  SCATTERING  SURFACE,  SATISFYING  THE  SCALAR  LEONTOVICH 
BOUNDARY  CONDITION  (ET  -  (n  •  ET)  n  )  =  p  Z  n  x  HT  . 

“  “  O  — 

Applying  scalar  and  vector  product  operations  of  equations  (4.3. 1)  to 

(4.3.4)  onto  one  another  a  complex  set  of  interdependent  scalar  and  vector 

equations  results  which  are  given  in  the  appendix.  Inspecting  these  equations 

it  was  found  practical  to  employ  the  following  three  pairs  of  vector  expressions 

T  T 

as  fundamental  vectors  for  further  analysis  (where  E  =  E  t  and  H  =  H  )  : 
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-i,  n 

=  A  +  A  =  [(EiEliiiflHiff]] 

(4.4.1) 

❖ 

■  -±I 

(4.4. la) 

-n 

■  "  -n 

(4.4. lb) 

Thus  Aj  and 

A  are  purely  imaginary  vector  quantities 

-I,n 

=  B  +  B  =  fj  (E*  x  H)  +  n*  (E  X  IP) 

)  (4.4.2) 

B* 

’  ’Si 

(4.4.2a) 

5‘n 

■ 

(4.4.2b) 

-i,  n 

=  C  +C  =  Q)  (£  X  H)  +  n'ME»x  H*j] 

(4.4.3) 

-i 

■  -  Si 

(4.4.3a) 

51 i 

5n 

(4.4.3b) 

In  addition  to  these  three  pairs  of  fundamental  vectors, 

introduce  the  following  notations: 

it  was  found  useful  to 

Bi 

*  AjXBj  =  D* 

(4.4.4) 

=  —  U  x  —  U  -  -  B* 

(4.4.5) 

E  =  A  xA  =  E  *  =  2AjXA2  =  2  r)  r)*(E  x  E*)  x  (H  x  H*)  (4.4.6) 

In  =  IIX1II  =  "^n=  2  B  x  B2  =  2  rj  rj* (E* x  H*)  x  (E  x  H*)  (4.4.7) 
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5m  ■ 

C  x  C  I;[  ■  -EJj  =  2  Cj  x  C2  =  2i)i)*(ExHlx(E*xH*l 

(4.4.8) 

l^j 

t-H 

II 

(4.4.9) 

F 

£II 

-  F*  =  BjXAjj 

(4.4. 10) 

Employing  properties  of  the  derived  set  of  scalar  and  vector  equations  as 
given  in  the  appendix,  three  orthogonal  vectors  out  of  the  thirteen  vectors 
defined  must  be  found  to  which  all  the  remaining  vectors  can  be  simply  associated. 
Inspecting  equations (4. 6. 6a,  9b,  13)  of  the  appendix,  two  of  the  required  vectors 
may  be  found  in  A ^  and  B  ,  since 

n  .  A  =  0  (4.4.11a) 

n  •  Bj  =  0  (4.4.11b) 


B i  =  n  x  A i  (4.4.11c) 

Aj  =  -nxBj  .  (4.4. lid) 


Thus  the  purely  imaginary  (4. 4.  la,  4.4.2a)  vector  quantities  A^  and  B^  are 
perpendicular  to  each  other  (4. 4. 11c,  4.  4.  lid)  and  in  fact  tangent  to  the  local 
scattering  surface  (4.4. 11a,  4.4.  lib)  .  Furthermore,  it  can  be  shown  that 


A  ^  .  A  j  =  (B  j  x  n)  •  A  j  =  (BjXn)  •  (BjXn) 


=  B  •  B  -  (n  •  B x )2  =  Bj  •  Bj 


(4.4. 12a) 


and 

A  ‘  B  =  0  .  (4.4.12b) 

Since  A^  and  B^  are  tangent  to  the  local  scattering  surface  and  perpendicular 
to  one  another,  its  cross  product  (4.4.4),  apurely  real  vector  quantity,  must 
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be  directed  along  the  local  outward  normal  n  of  the  scattering  surface,  where 


2i  ■  A,  -  Bj  ■ 

-  (Aj  •  A^  n  =  -  (Bj 

• B^ n  = [ 

_S  •  (A  j  x  Bj)] 

n  (4.4.13) 

and  with  equation  ( 

.  .  ),  where  it  was  shown  that: 

[n  •  (AjXBj)] 

■  - 2  V  c2  *  -2(4, 

•  a2)  +  b,- 

V 

(4.4. 14a) 

-I  -  -2^r  C2 

-  -2i  (Ar^  +  sr 

V 

• 

(4.4.14b) 

Since  no  further  triplet  of  orthogonal  vectors  was  found  and  the  remaining  vector 
quantities  as  defined  in  (4 . 4. 1)  to  (4.4. 10)  can  be  uniquely  decomposed  along 
the  directions  of  A  ,  Bj,  and  D  ^  ,  the  following  theorem  can  be  formulated: 

Theorem:  If  the  electromagnetic  behavior  in  the  vicinity  of  a  scatterer  satisfies 
the  Leontovich  Boundary  condition 

ET  -  (n  *  ET)  n  =  q  Z  nx  HT 

—  o  — 

then  the  following  two  purely  imaginary  vectors 

T  T*  T  T* 

A  x  =  (E  x  E  )  -  rj  rj*(H  x  H  ) 

and 

T*  T  T  T* 

B  =  n  (E  x  H  )  -  rj*  (E  x  H  ) 

are  orthogonal  vector  quantities  which  lie  in  the  local  tangent  plane  of  the  scatterer 
and  its  cross  product,  a  purely  real  vector  quantity, 

=  4jx®j  =  n  [q  •  ^AjXBj)] 

is  directed  along  the  outward  local  normal  of  the  scatterer.  The  absolute  values 
of  these  three  vector  quantities  are  identical,  where 

Aj  •  Aj  =  Bj  •  Bj  =  -  (n  ■  Dj) 
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and  Aj  ,  B  ,  Dj  constitute  a  right-handed  orthogonal  vector  triplet. 

These  peculiar  properties  result  in  the  following  characteristic  equations 
from  which  a  set  of  inverse  scattering  boundary  conditions  may  possibly  derived. 

Since  with  (4.4.13) 

=  Sj  Sj  =  -(S-Dj)  =  +8^.^ 

=  2<W  W 

two  characteristic  equations  for  the  determination  of  the  amplitude  and  the  phase 
of  the  relative  impedance  rj  are  found: 

(Ex  E*)  •  (E  x  E*)  +  (r/  rf*)2  (H  x  H*)  .(Hxlf)  =  -  2  r?  if(E*xH)  •  (E  x  H*) 

(4.4.15a) 

and 

rj2  (E*x  H)  •  (E*x  H)  +  rf &  (E  x  H*)  •  (E  x  H*)  =  -  2r\  rj*( E  x  E*)  •  (H  x  H*).  (4.4. 15b) 

For  the  derivation  of  (4. 4. 15a)  and  (  4. 4. 15b),  equation  (4.4. 12b)  was  not  used, 
thus  the  orthogonality  condition  of  A  and  B  may  be  employed  as  the  character¬ 
istic  equation  for  the  determination  of  the  surface  locus,  i.e. 

A  '  B  =  0  or  (4.4.16a) 

[(E  x  E*)  -  n  rfHH  x  H*)J  •  [rj  (E*x  H)  -  rj*( E  x  H*)]  =  0 
which  can  be  rewritten  into 

[(E  •  E*)  -  rj  rj*( H  •  H*)]  [rj  •  (E*  •  H)  +  tj*(E  •  H*)]  = 

=  n  (E  •  H)  [(E*  •  E*)  -  rf2  (H*  •  H*)]  +  rj*(E*  •  H>:<)  [(E  •  E )  -  r]2  (H  •  H )]  . 

(4.4.16b) 

Although  (4.4. 16)  has  not  been  employed  in  the  derivation  of  (4.  4. 15a,  b),  these 
three  equations  may  be  linearly  dependent  which  still  must  be  investigated. 

The  fourth  characteristic  equation  for  the  determination  of  the  surface 
normal,  then  is  given  by 
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4 

n 


(4.4.  17b) 


=  T 


Bl 


A,  x  B 


n  *  —  I 


6  '  [^lx5:] 


ygt 

2  Si  •  e2 


-  |  £(E  x  E4)  -  n  n”(H  x  H4)J  x  <E4x  H)  -  r74(E  x  H4)J^ 
2  n  iftti  H)  •  <E*x  H4) 


which  is  dependent  upon  rj  and  so  far  could  not  be  formulated  solely  in  terms 
of  the  field  quantities. 

Inspecting  these  characteristic  equations  it  is  found  that  a  fourfold 

solution  for  the  relative  impedance  rj  results  as  was  predicted,  furthermore 

these  rj  represent  a  set  of  four  complementary  solutions.  For  the  degenerate 

case  r)->0,  equation  (4.4. 15a)  degenerates  into  the  inverse  scattering  boundary 
X  T* 

condition  E  x  E  =0  which  was  derived  by  Weston  in  7644-  1-F  (chapter  X). 

To  determine  the  proper  amplitude  and  phase  of  r)  ,  the 
formulation  (4.3.1)  of  the  Leontovich  equation  must  be  satisfied  which  would 
require  a  characteristic  equation  for  the  surface  normal  n  being  independent 
from  T]  ,  otherwise  a  unique  solution  may  not  be  found  from  the  orthogonality 
condition  (4.4. 16a).  Namely  were  r?  a  function  of  frequency  then  also  n  may 
become  frequency  dependent.  However  it  may  be  possible  to  express  D ^  in 
terms  of  cross  procuts  of  A  9,  B  „  or  C.  „  such  that  the  normal  becomes 
a  closed  form  expression  in  terms  of  the  field  quantities  only. 

Although  the  characteristic  equations  (4.4. 15,  16  and  17)  cannot  be  applied 
in  the  general  problem  of  inverse  scattering  for  the  determination  of  the  constituent 
parameters  S  (R,  B,  0),  n  (R,  6, 0)  and  n  (R,  6,  0)  due  to  the  inherent  vector  nature, 
the  presented  derivation  suggests  that  such  a  formulation  of  local  inverse  scattering 
boundary  conditions  may  exist.  A  more  sophisticated  analysis,  employing  the 
vector  form  of  the  Leontovich  boundary  condition,  may  yield  the  desired  character¬ 
istic  equations. 
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4.5  SUMMARY 

It  was  shown  that  it  is  possible  to  derive  a  set  of  scalar  and  vector  equations 
from  the  complementary  set  of  the  four  formulations  of  the  Leontovich  boundary 
condition  which  may  be  employed  as  local  boundary  conditions  to  the  problem 
of  continuous  wave  inverse  scattering.  Although  the  derivations  are  based  upon 
the  formulation  of  a  scalar  Leontovich  impedance,  this  first  approach  reveals 
a  rather  important  result.  Namely  inspecting  the  characteristic  equation  (4.4. 15a) 
it  can  be  shown  that  the  associated  scattering  surface  locus  can  be  found  if  solely 
the  absolute  value  of  the  impedance  is  known.  Thus  if  the  material  character 
of  the  surface  is  known,  the  following  local  inverse  scattering  boundary  conditions 
may  be  applied  to  determine  the  surface  locus  S(R,0,  0),  i.e.  the  shape  of  the 
scatterer: 


a)  ri  =  0  (perfectly  conducting  scatterer): 

T  T*  T  T * 

(E  x  E  )  •  (E  x  E  )  =  0  (4.5.1) 

0)  r\  =  oo  (perfect  magnetic  conductor  or  scatterer  with  plasma  coating): 


T  T*  T  T * 

(H  x  H  )  •  (H  x  H  )  =  0 


(4.5.2) 


7)  0  <  r)  <  00,  /r)  |  /  1  (conducting  scatterer  for  which  only  the  absolute 

value  of  the  relative  Leontovich  impedance  has 
to  be  known): 

(ET  x  ET*)  •  (ET  x  ET*)  +  (nrf)2  (HT  x  HT1  *  (HT  x  HT*)  = 

rr#  t  rp  rr# 

=  -  2  rjr)*  (E 1  x  H  )  •  (E  x  H  )  .  (4.5.3) 


These  boundary  conditions  are  necessary  but  not  sufficient, however  they  are  local 
conditions  which  permit  a  point  by  point  determination  of  the  surface  locus  and 
its  application  of  two  or  more  different  operational  frequencies  may  result  in 
the  unique  specification  of  the  proper  locus  out  of  the  infinite  set  of  the  obtained 
loci  as  was  discussed  in  detail  in  section  (4.2). 
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Furthermore  the  properties  of  the  formulated  theorem  of  section  (4.4) 
suggest  the  formulation  of  a  more  general  set  of  local  inverse  scattering  boundary 
conditions  for  which  the  a  priori  knowledge  of  the  relative  Leontovich  impedance 
is  not  required.  Such  a  formulation  requires  first  of  all  the  proof  that  the  charac¬ 
teristic  equations  (4.4. 15a),  (4.4. 15b)  and  (4.4. 16a)  are  linearly  independent. 
Then  a  formulation  of  the  surface  normal  must  be  sought  which  is  independent 
of  r)  so  that  the  impedance  corresponding  to  the  proper  plane  can  be  selected. 

So  far  no  explicite  physical  interpretation  on  the  nature  of  these  boundary 
conditions  was  presented,  preliminary  investigations  however  suggest  that  a 
relationship  between  the  electromagnetic  stress  tensor  (J.  A.  Stratton,  1941) 
and  the  properties  of  the  above  state  theorem  exists.  It  is  therefore  proposed 
to  investigate  such  existing  relationships  in  all  detail. 

4.6  APPENDIX 

T 

For  convenience,  the  total  magnetic  field  vector  H  =  H  is  normalized 

T 

with  respect  to  the  total  electric  field  vector  E  =  E  by  the  intrinsic  impedance 
of  free  space  Zo  =  120  7 ,  where  Eqs.  (4.3. 1)  to  (4.3.4)  will  assume  the 
following  form 


■J_E-  (n  •  E)  n] 

/  * 

-  r)  (n 

x  H)j  = 

0 

I 

|(Exi\)-ii  [h  ■ 

-<n  •  H)  n]  j  = 

0 

II 

f  [e*  -  (n  •  E*)  n] 

|  -i)*  (nx  H*)j 

!>  =  o 

in 

^  ( E  *  x  n )  -  r]  *  1 

-  ( n 

•  H*)  n 

]]  =0. 

IV 

The  scalar  and  vector  product  operations  of  these  equations  onto  another  result 
into  the  following  relationships: 

4.6.1  SCALAR  PRODUCTS 

I  •  I  or  II  •  II 

(E  •  E)  -  n2  (H  •  H)  =  (n  •  E)2  -  p2  (n  •  H)2  (4.6.  la) 
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(4.6.1b) 


rj  [e  •  (n  x  H)]  =  (E  •  E)  -  (S  •  E)2  =  n  [h  ’  (E  X  n)]  = 

2  (H  •  H)  -  n2  (n  •  H)2 
or 

2r7[S-(HxE)]  =  [(E  •  E)+ n2(H  •  H)]  -  [(n  •  E)2+ n2(n  •  H)2]  (4.6.1c) 

i*  •  i*  or  11*  •  n* 

(E*  •  E*)  -  rf£  (H*  •  H*)  =  (n  •  E*)2  -  n*2  (n  •  H*)2  (4.6.2a) 

rj*  [|n  •  (H^xE*)]  =  rj*  E*  •  (n  x  H*)J  =  n *  H*  •  (E*xn)J  = 

=  (E*  •  E*)  -  (n  •  E*)2  -  r)  '2  ( H*  •  H*)  -  n*2  (a  ■  H*)2  (4.6.2b) 

or 

2 n*  [n  •  (H* x  E’^]  =  [(E*  •  E*)  +  n*2  (H*  •  H*)]  -  [(n  •  E*)2  +  rf *2  (n  •  H>lc)2] 

(4.6.2c) 

I  •  n  and  I*  -  II* 

rj(E  •  H)  =  rj(n  •  E)  (n  •  H)  (4.6.3a) 

17* (E*  •  H*)  =  n>;t(n-  E*)  (n  •  H*)  (4.6.3b) 

I  .  I* 

(nx  E)  •  (nx  E*)  -  r)rj*(nx  H)  •  (n  x  H*)  =  0  (4.6.4a) 

or 

(E  •  E*)  -  nnHll  •  H*)  =  (n  •  E)  (n  •  E*)  -  nrfH n  •  H)  (n  •  H*)  (4.6.4b) 

II  •  n* 

n-  frj(E*x  H)  -  rj*(E  x  H*)  }  =0  (4.6.5) 
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<  cl 


I  •  IF 


•  £(E  x  E*)  -  r)r?*(Hx  H*)]  =  0 
rj  (n  x  E*)  •  (nxHj  +  ij^nxE)  *  (n  x  H*)  =  0 


or 


T)  (E*  •  H)  +  n*( E  -IF}  =  n( n  •  E*)  +  77*  ( n  •  E)  (n  •  H*) 

4.6.2  VECTOR  PRODUCTS 
I  x  I  and  II  x  II 


(4.6.6a) 

(4.6.6b) 

(4.6.6c) 


2 

(E  x  E)  =  r)  (HxH)  =  0 

(4.6.7a) 

F  x  F  and  IF1'  x  II* 

2 

(E*xE*)  =  n*  (H*x  IF)  =  0 

(4.6.7b) 

I  x  II  and  F*  x  IF15 

n  (E  x  H)  =  ^[r?  n  •  (H  x  E)]  -  [(E  •  E)  +  r?  (H  •  H)]j  n  +  [(n  • 

E)  E  +  n2  (n  •  H)  h] 

=  [n2  (H  •  H)  -  (n  •  E)2]  n  +  [(n  •  E)  E  +  r}2  (n  •  H)  H 

]  = 

=  [(E  •  E)  -  n2  (n  •  H)2]  n  +  [(n  *  E)  E  +  r?  (n  *  H)  h] 

(4.6.8a) 

n*(E*x  H*)  =  (|>n  •  ( H* x  E*)]  -  [(E*  •  E*)  +  rff2  (H*  •  H*)]] 

A 

n  + 

+  f(n  •  E*)  E*  +  n*2  (n  •  H*)  H*1  = 

=  [rr  (H*  •  H*)  -  (n  •  E*)2  ]  n  +  [(n  •  E*)  E*+  n*(n  •  H*)  H*]  = 

=  [(E*  *  E*)  -  n*2  (n  •  H*)2  J  n  +  [(n  •  E*)  E*+  r?*2  (n  •  H':<)  H*^] 


(4.6.8b) 
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IxF  and  II  x  II* 


n  •  [(E  x  E*)  +  r)T]*( H  x  H*)^  =  [rf‘(n  x  E)  •  (n  x  H*)  -  r)  (n  x  E*)  •  (n  x  H)] 
since 


(4.6.9a) 


rj  (n  x  E*)  •  (n  x  H) 


n  •  (E  x  E*)  =  -  r)r7*n  •  (H  x  H*) 


(4.6.9b) 


r)*(n  x  E)  #(n  x  H)  =  n  •  (E  x  E*)  =  r)rffn  •  (H  x  H*)  (4.6.9c) 

thus 

n  •  [(Ex  E*)  -  r/rr(H  x  H*)]  =  0  (4.6.9d) 

and 

n  [rj  (n  x  E*)  *  (n  x  H)  +  r]*(n  x  E)  •  (n  x  H*)"]  =  0  (4.  6. 9e) 

I  x  II* 

n  [(n  x  E)  •  (n  x  E*)  -  ryri'H n  x  H)  •  (n  x  H*)]  =  0  (4.6. 10a) 

n  |(nxE)  •  (nx  E*)  +  r)  r)*  (n  x  H)  *  (n  x  H*)|  = 

=  rj*  n  x  [(E  x  H*)  x  n  [  +  r)  (n  x  H)  x  (n  x  E*)  -  rj* ( E  x  H*)  (4.6. 10b) 

II  x  I* 

n  |(n  x  E)  •  (n  x  E*)  +  r)rj*(n  x  H)  *  (n  x  H*)  j  = 

=  r)*(n  x  H*)  x  (n  x  E)  -  r)  (E*x  H)  +  t)  n  x  [ ( E * x  H)  x  n  ~\  (4.6. 10c) 

(I  x  IT*)  -  (II  x  P'<) 


[rj  (E*  x  H)  -  rj*  (E  x  H* ) ~j  =  r)*(n  x  H*)  x  (n  x  E)  -  rj  (n  x  H)  x  (n  x  E*) 

+  rj  n  x  [(E*x  H)  x  n[J  -  77*  n  x  [(E  x  H*)  x  n  "j  = 

=  nx  [(Ex  E*)  -  r)r)*(H  x  H*)J  .  (4.6.11) 
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Equation  (4.6.11)  constitutes  an  important  result  namely  that  a  set  of  orthogonal 
vectors  is  found  which  are  tangent  to  the  local  scattering  surface  since  by 
(4.6.5)  and  (4.6.6a) 

n  •  [rj  (E*x  H)  -  n*( E  x  H*)  ]  =  0 

and 

n  •  [(Ex  E*)  -  rjrj*(  H  x  H*)]  =  0 


V erification  of  ( 4 . 6 .  1 1 ) 

[n  ( E * x  H)  -  rJ*(E  x  H*)]  =  [e*x  (E  X  n)  -  E  X  (E*xn)] 

+  rjq*  ^(n  •  H)  [(n  x  H*)  x  n  J  -  (n  •  H*)  [(n  x  H)  x  n[  j  = 

=  [(E*  .  n)  E  -  n  (E  •  E*)  -  (E  •  n)  E*  +  n  (E  •  E*)] 

+  nr)*  [(H  •  n)  H*  -  (n  .  H)  (n  •  H*)  n  -  (n  .  H*)  +  (n  .  H*)  (n  •  H)  n] 

=  [(E*x  E)  x  n_j  +  r)rj*  [(H  x  H*)  x  n  [  - 
=  n  x  [(E  x  E*)  -  r;r)*(H  x  H*)J 
There  exists  another  formulation  of  (4.6. 11)  given  by 

[(E  x  E*)  -  rjn*(H  x  If*)]  =  -n  x  [r?  (E*x  H)  -  n*(E  x  H*)[|  (4.6.12) 


which  may  be  verified  as  follows: 

[(E  x  E*)  -  r)if;‘(H  x  H*)[  =  [e  x  (n  x  H*)  j  -  q  [hx(E,jx  n)"]^ 

+  |(n  .  E*)  [*] E  x  n^  -  r)rj*(n  •  H*)  [h  x  n]  j  = 

=  n)  [j*(E  •  H*)  +  rj  (E*  •  H)]  -  [n*(n  •  E)(n  •  H*)  +  rj  (n  •  E*)  (n  •  H)]  f 
v - v - ' 

=  0,  according  to  (4. 6. 6c) 

+  |[n  (n  •  E*)  H  -  n  (n  •  H)  E*]  +  [rj*( n  .  H*)  E  -  rj*(E  •  n)  H*]  ’  = 

=  r\  [n  x  (H  x  E*)j  +r)*[nx(Ex  H*)"]  = 

=  -  n  x  [rj  (E*x  H)  -  rj*( E  x  H*)  J  q.e.d. 
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<  a 


This  concludes  the  derivation  of  scalar  and  vector  product  operations  which 
resulted  into  the  following  identities: 

{[<E  x  E*)  -  n T7* (H  x  H*)]  +  n  x  [rj  (E* x  H)  -  rf( E  x  H*)]  ]  =  0  (4.6.  13a) 

or 

/ 

'  [(E  x  E*)  -  nn*(H  X  H*)]  X  n  +  [rj  (E*x  H)  -  n*(E  x  H*)]j  =0  .  (4.6.13b) 

4.6.3  DERIVATION  OF  ADDITIONAL  VECTOR  IDENTITIES 

It  can  be  shown  that  D ^  =  n  (Dj  •  n)  can  be  reformulated  as: 

•  D:  -  n  •  (A  x  B  )  =  -  2  |  [nn*( n  •  E)  (n  •  E*)  (H  •  H*) 

+  rin’Hn  •  H)  (n  •  H*)  (E  •  E*)]  -  [riri*(n  •  H)  (n  •  E*)  (E  •  H*) 

+  nn*( n  •  H*)  (E  •  n)  (E*  •  H)])  = 

=  -  2  nn*  [(E  X  H)  •  (E*x  H*)  ]  =  -  2  r\rf  [(E  x  E*)  •  (H  x  H*) 

+  (E*x  H)  •  (E  x  H*)j  =  -2  nn*  [(E  •  E*)  (H  •  H*)  -  r)rf( E*  •  H)  (E  •  H*)] 

(4.6.14) 
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V 


MATRIX  INVERSION  FOR  THE  CASE  OF  END-ON 
INCIDENCE  ONTO  A  ROTATIONALLY  SYMMETRIC  SCATTERER 

5.1  FORMULATION  OF  THE  MATRIX 

In  problems  of  C.  W.  inverse  scattering  it  is  assumed  that  for  a  given 
incident  field  the  measured  far  scattered  field  can  be  obtained  for  a  sufficiently 
large  number  of  bistatic  angles.  If  the  direction  of  the  incident  wave  is  chosen 
along  and  in  direction  of  the  negative  z  -  axis  of  a  spherical  coordinate  system 
and  the  polarization  vector  e  of  the  transmitted  wave  is  parallel  and  in  direction 
of  the  positive  x  -  axis,  then  the  scattered  field  may  be  represented  by  a  series 
expansion  into  vector  spherical  wave  functions: 


oo 


n 


(5.1.1) 


n=l  m=0 


where  the  vector  spherical  wave  functions  are  defined  by: 


A 


e 


(5.1.2) 


e 


A 


R  h  (kR) 
n 


0  (5.1.3) 


A 


with 


+  m  sin  6  Pm(  cos  6) 
n 


(5.1.4) 
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Rm(0)  = 
n 


m, 

3  P  (cos  0) 
n _ 

30 


=  ~  ["(n-m+l)  (n+m)  Pm  ^cos  0)  -  pm+1(cos  0) 

2  L  n  n  J 


.(5.1.5) 


(5.1.5) 

For  a  perfectly  conducting  sphere  the  expansion  coefficients  aemm  l 

o 

be^n  bear  the  following  relationship  with  the  expansion  coefficients  given  by 
o 

Stratton  (1941),  where  m  =  1  for  nose-on  incidence  onto  a  rotationally  symmet¬ 
ric  body: 


a 


oln 


i(-l) 


n+ 1  (2n+  1) 
( h  +  l)n 


a 

n 


(5.1.6a) 


b 


eln 


i  (-D 


(2n+  1) 

( n+  1)  n 


b 

n 


(5. 1.6b) 


Employing  the  asymptotic  approximation  of  the  Hankel  functions 


lim 
x  ->  oo 


,(1)(x)l  =  (-l)1**1  SSEi* 
n  x 


(5. 1.7a) 


lim 
x  ->  oo 


_d_ 

dx 


[x  ^U(x)] 


=  (-i) 


n 


expix 

x 


(5.1.7b) 


and  extracting  the  factor  [  eX^ 1  x  j  from  the  measured  field  components,  the 
far  scattered  field  is  given  by 

N 

Ee  0C>  =  cos  0 

c  n=l 

N 


[■ 


a,  S1(cos0)  +  b1  R^cos  0  )  1  (5.1.8a) 

oln  n  c  eln  n  c  J 


e!  (e  J  ) 

0  c  r  c 


R^(cos  0  )  +  b  , 
n  c  eln 


S1  (cos  0  )  I 
n  c  J 


(5.1.8b) 


where  N  denotes  the  total  number  of  receiver  aspect  angles,  and  for  the  c 
receiver  the  field  components  are  given  by 
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ES  (0  ,  0  )  ,  E®  (0  ,  0  )  . 

0  c  c  0  c  c 
c  c 

For  the  formulation  of  a  suitable  near  field  representation,  the  unknown 

expansion  coefficients  a  ,  b  must  be  computed.  To  do  so.  the  field 

oln  eln  .  . 

s  s  11 

components  ,  E ^  ,  and  the  inherent  terms  S  (0),  R  (0),  cos  0  must 
0  0  n  n  rc 

c  c 

be  arranged  in  a  matrix  formulation  so  that  the  most  suitable  matrix  repre¬ 
sentation  results  which  satisfies  stability  criterions,  encountered  in  the  inver¬ 
sion  of  such  a  matrix. 

Inspecting  Eqs.  (5. 1.8a,  b)  which  represent  the  basis  of  further  analysis, 
these  equations  may  be  rewritten  into  the  formulation 


E  =  A  X 


(5.1.9) 


T 

where  the  transpose  E  of  E  is  given  by: 


E9.  (VV-  V92'  Ve3'*3>‘  •  •  • 


91 


•  •  •  E»lWE(i1<,i^i)EAV'  ■ 

N  1  2 


E j  (0  .  0  ) 

0  '  N  PN 
rN  -1 


(5.1.10a) 


which  consists  of  2N  complex  elements,  and  so  does  the  transpose  X  of 
X  which  represents  the  unknown  coefficients: 


x  :  <a  b„ ii>  a  io«  b„io*  ■  •  •  a  b  )  .  (5.1.10b) 

oil  ell  ol2  el2  oln  eln 

With  this  arrangement  the  matrix  elements  of  A  are  defined  as  well,  where 
th 

the  r  row  for  1  <  r  <  N  is  given  by 

S \  (cos  0  )  cos  0  ,  R*  (cos  0  )  cos  0  ,  (cos  0  )  cos  0  ,  R^  (cos  0  )  cos  0  , 

LI  r  r  1  r  r  2  r  rr  2  r  rr 

...  S*  (cos  0  )  cos  0  ,  R*  (cos  0  )  cos  0  1  (5.1.  11a) 

N  r  r  N  r  rrJ 


42 


and  the  (N  +  r;  row  for  1  <  r  <  N  is  given  by 

f-  R^(cos  0  )  sin  0  -  S^(cos  0  )  sin  0  ,  -  R^(cos  0  )  sin  0  , 

LI  r  rr^  1  r  rr  2  r  rv 

(5.1. lib) 

-  S*  (cos  0  )  sin  0  ,  .  .  .  ,  -  R* (cos  0  )  sin  0  ,  -  S*(cos  0  )  sin  0  "1 
2  r  r  N  r  r  N  r  rj 

The  resulting  matrix  A  can  however  be  decomposed  into  the  product  of  two 
square  matrices: 

A  (0  ,  0  ,  n)  =  B(0  )  C(0  ,  n)  .  (5.1.12) 

c  c  c  c 

The  properties  of  this  matrix  must  be  investigated  in  detail  to  gain  more  informa¬ 
tion  about  stability  criterions . 

5.2  PROPERTIES  OF  MATRIX  A  (0  ,  0  ,  n)  =  B(0  )  C(0  ,  n) 

c  c  c  c 

Matrix  A  is  a  square  matrix  with  2N  x  2N  purely  real  elements  for 

real  angles  0  ,  and  0  .  Matrix  B(0  )  is  diagonal,  solely  incorporating  the 

c  c  c 

0  -  dependence  and  its  elements  are  given  by: 


b  =  cos  0 
rr  rc 


for  1  <  r  <  N 
'  c 

^H-r,  N+r 

The  determinant  of  B  is  always  smaller  than  unity,  where 


sin  0  for  N  <  N+r  <  2N 
c  — 


(5.2.1a) 

(5.2.1b) 


Det 


fB) 


(-1) 


N  N 


.N 


TT  sin  2  0C 

c  -  1 


(5.2.2) 


This  matrix  will  become  singular  for 


0  =  0  +  m  ir/2  ,  m  =  0,  1,  2 

c 


(5.2.3a) 


and  is  maximum  for 


0^  =  7r/4  +  m  7r/2 


(5.2.3b) 
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and  it  can  be  seen  that  otherwise  no  restrictions  are  imposed,  e.g.  the  choice 
of  0  =  (5  :  fl  =  .  ,  .  :  f)  is  permissible.  However,  it  may  be  advisable 

C  1  ^  IN 

to  extract  the  0  -  dependence  entirely  from  the  field  components  before  in¬ 
verting  the  matrix,  since  the  0  -  dependence  solely  reduces  the  magnitude 
of  determinant  Det  |A  j  . 

Matrix  C  (0  )  ,  solely  incorporating  the  0  -  dependence  in  terms  of 
c 

Legendre  functions,  needs  further  thorough  investigation.  Inspecting  Eqs. 

(5.1. 8a,  b)  it  is  recognized  that  always  two  of  the  2N  x  2N  elements  are  identi¬ 
cal,  and  for  even  N  matrix  C  can  be  decomposed  into  the  sum  of  two  matrices 
each  containing  non-identical  coefficients,  where  for  N  =  2  we  find: 


C  = 


C11 

C12 

C13 

C14 

C11  ° 

C13 

0 

C21 

C22 

C23 

C24 

=  C  +  C„  = 

°  C22 

0 

C24 

C31 

C32 

C33 

C34 

1  2 

C31  ° 

C33 

0 

C41 

C42 

C43 

C44 

°  C42 

0 

1 

O 

0  °12  °  c14 

C21  °  °23  0 

0  °32  0  °34 

°41  °  C43  ° 


(5.2.4) 


where 


32  =  S!<V 

=  1 

C13  =  c34 

:42-Si(e2> 

=  1 

C23  =  C44 

3i=Ki<v 

=  cos  0^ 

C14  =  C33 

41=R!(92) 

=  cos  0 

o 

to 

II 

o 

00 

S^)  -  3  cos  e2 
R*(0  )  -  3  cos  2  0 

Ct  L  1 

rJ(02)  =  3  cos  2  0, 


(5.2.5) 


Furthermore,  it  can  be  shown  that 


44 


Detfr]  "  DetlC2i 


(5.2.6a) 


since 


T  C  R 


(5.2.6b) 


where 


T  = 


’0 

1 

0 

0" 

’0 

0 

1 

(T 

1 

0 

0 

0 

R  = 

0 

0 

0 

1 

0 

0 

0 

1 

1 

0 

0 

0 

_0 

0 

1 

0_ 

_0 

1 

0 

0_ 

(5.2.6c) 


This  property  may  be  employed  to  derive  the  determinant  of  C  more  easily, 
where 

C=C1  +  C2  =  CjU+C"1^)  -  (C1C21  +  I)  C2  (5.2.7a) 


where 


Det  {C  j  =  Det  ^  J  Det  ^1  +  C~ 1  C2  } 


(5.2.7b) 


For  the  chosen  case  N  =  2,  the  determinants  result  into: 

““‘K  }  -  (c12  °13  -  'll  c14)  (°22  °23  ‘  °21  °24) 


-  (3  sin^0  )  (3  sin^09  ) 


Det  |l  +  c”1  C2  |  =  [2  (cos  @1  -  cos  02)  J 


(5.2.8a) 

(5.2.  8b) 


thus 


Det  |c  j  =  -  [j3  sin^Q^  ][  3  sin^02j  £2  (cos  0^  -  cos  @2)j  -(5.2. 

N-2 


8c) 
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Employing  this  scheme,  the  determinants  for  the  cases  N  =  4  and  N  =  6  were 
evaluated.  There  it  was  found  convenient  to  formulate  the  determinants  into  a 
continued  fraction  expansion  for  increasing  N,  where 


Det 


ft- 


C1  sin  0 


1 

c2sin202 

jjcos  0^  -  cos  0  J 

- 1 

1 

c38to2e3| 

COS  0g  -  COS  02^{c°S  0g  -  cos  0^ 

2 

C  .  sin  0 A 
4  4 


cos  0 ^  -  cos  0gj^jcos  0 ^  -  cos  0  2]2[o°s  0^  -  cos  0 


J 


CL  sin  0 


5  5L  5 


cos  0„  -  cos  0  J  [cos  -  cos  0  J  [cos  0^  -  cos  |cos  0^  -  cos  0  J " 


CL  sin  0, 
b  ( 


cos  6„  -  cos  0_ 
b  aJ 


j^COS  0g  -  cos  0 


etc. 


CN  Sl"  % 


cos  0XT 

-  COS  0VT 

L  N 

N  -  1 

u 

L  N  lj 

N 


' 


(5. 2.9a) 
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The  constant  multipliers  were  found  to  be  : 


c  =  i  =  (-i)1  +  1  -  •  i2 
1  1 


C2  =  -22  32  M-1)2+1^/  22 


e3.^- (-D3+  1  ‘2  f  f  -  32 
3  2  1  •  2^ 


C  4---22.S2-72 I2//-/-*2 
4  1  •  2  •  3 


_  3  •  54-  72  5+  1  1  •  32-  52-  72  -  92  2 

c - g - (-1)  •  - 2 - 2 - 7 - 2 - 5 

D  2  1  •  2  •  3  •  4 


622  „  2  2  2  2  2  2 
„  _  -3  •  7  •  11  _  ,  lX6+  1  1  •  3  •  5  •  7  •  9  •  11  c2 

e6  ”  4  “  '  2  2  2  2  2  ’  b 

2  1  •  2  •  3  •  4  •  5 


C  =  (-l) 
c 


c  +  1  [(2c  -  1)  '■  l] 

[(c-l)O  2 


(5.2.9b) 


(2c -1)!'.  =  1  •  3  •  5  -  7  -  9 - (2c -1) 

where  computational  results  are  given  in  Table  V-l.  Inspecting  (5.  2.  9)  it  can 
be  seen  that  Det  jc  ( 0 2N  x  2N)J  can  be  decomposed  into  the  product  of  two 
N  x  N  determinants: 


C  (0  j  2  N  x  2N)  =  D  (sin20  ,  Nx  N)  •  E2  (cos  6  ,  Nx  N) 
c  c  c 

D  is  a  diagonal  matrix  with  diagonal  element  of  the  form 

fr20-.1*'-'-]2 


(5.2.10a) 


d  =(-l) 
cc 


c+  1 


[(c-1)'.] 


9  c2sin20 
2  c 


(5.2.10b) 
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and  E  is  of  the  Vandermonde  type,  where  Vandermonde's  determinant  is 
here  given  by: 


V=E  = 


2 

lcos  ©^cos  0  •  • 

N  -  1 
•  ’  cos  0^ 

( 

,  n  2n 

lcos 0  cos  0 

N-  1 

cos  0 

Ci 

\ 

2 

N  - 1 

lcos 0  cos  0 

cos  0 

1 

N-l 

/ 

1COS0 . 

cos  0 

! 

N 

N 

i  V 

(cos  0 -cos  0  )(cos  0  -  cos  0 )•  •  •  (cos  0  -cos  0  ) 

u  1  j  1  IN  1 


(cos  0 -COS  0  )(cOS  0  -COS  0  )•  ••  (cos  0  -cos  0  )  I 
o  2  4  2  N  2 


(cos  0  —  cos  0  ,) 

N  N  -  1 


N 


N 


JJ  (cos0  -cos0)=  7J  2  sin  (  -S--  —  )  sin  (  r)  .(5.2.10c) 
N>r >  s> 1  r  S  N>r >  s> 1 


0+0  0-0 
3 _ 

2 


Thus  instead  of  the  continuous  fraction  expansion  a  closed -form  solution  for 
Det|c(0c,  N)^  can  be  given,  where 


Det 


c(0 , n)|  =  TT  (-dc+1 

C  ;  c  =  l 


c  (2c  -  1)'. 

—7 - tt: -  Sin  0 

(c-1)'.  c 


2  N 

n(cos 0  -cos  0  7 

T*  c 

N>r>  s> 1 


N  „  N 

-  fcotasj  rr 

c  =  l  N>r>  s  > 1 L 


•  0+0  0-0-12 
2(2s+l)  .  ,  s  rx  ,  ,  s  rx 

-j-  sm(— )sW— ) 


N 


(5.2. 11) 


Inspecting  (5.  2. 11)  it  is  found  that  for  a  set  of  N  aspect  angles  R  =  2  c  roots 

2  C-1 

will  be  encountered,  where  N  of  these  roots  are  the  sin  0  terms  and  the 

c 

remainder  is  the  product  of  the  squares  of  the  differences  of  the  argument  cos  0^ 

of  the  Legendre  functions,  for  all  mutations  of  the  N  polar  angles  0  .  Thus  the 

following  restrictions  must  be  imposed  onto  the  0  -dependence: 

c 
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i)  8  41  0  or  7r 
c 


ii)  8  +  6  ,  ,  c  =  2,  3,  .  .  .  .  N 

c  c-1 


iii)  The  distribution  of  the  polar  angles  8 must  be  chosen  so  that  the 


determinant  neither  becomes  too  much  nor  too  large  but  of  the  order  of  unity 
where  for  each  different  number  N  an  individual  optimization  procedure  must 
be  employed. 

5.3  SUMMARY 

The  obtained  results  may  be  summarized  in  the  following  theorem: 
Theorem:  V-l 

If  an  expansion  of  the  scattered  field  into  spherical  vector  wave  functions 
is  employed  for  the  first  dipole  case  (m  =1),  i.  e. ,  end-on  incidence  onto  a 
rotationally  symmetric  and  perfectly  conducting  scatterer,  then  the  determinant 


Det 


/a  (0,0;  N) 
c  c 


M  N 

=  (-nN  TT  (-i) 


yc  c  •  (2c-  l)1.1. 

2  JL  (c-1)’ 


2 

TT  (cos  6  -cos  8  ) 

'  '  ,  r  s 


c  =  1 


of  the  associated  scattering  matrix  ^(0^,  become  singular  for 


i)0=P77*  p  =  0  +  1,  1  2,  t  3 

c  ^ 


ii)0=q7r  ,  q  =  0  ,  1 


c 


iii)  8=8^.  ,  0  =  1,  2,  ...  (N  -  1) 

c  c+ 1 
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In  addition,  pseudo-singular  behavior  may  be  encountered  even  though 
(i),  (ii)  and  (iii)  are  not  satisfied  if  the  computational  aspect  angles  exclusively 
lie  within  narrow  cones  about  the  z-axis  of  the  computational  coordinate  system, 
thus  reducing  the  value  of  the  determinant  to  almost  zero. 

These  properties  are  illustrated  in  Figs.  5-1  and  5-2  for  transmitter- 
receiver  configurations  which  may  occur  in  practice  most  frequently. 

For  both  cases  the  receiver  aspect  angles  are  assumed  to  be  well 
distributed  within  a  narrow  cone  whose  invariant  axis  a  is  oriented  perpen¬ 
dicular  to  back  scatter  direction.  For  simplicity  it  is  assumed  that  0c  =  7t/4 

or  0  =  7t/4  +  7T,  and  main  attention  will  be  attributed  to  the  dependence  of  the 
c 

polar  angle  0  . 

0 

The  determinant  associated  with  the  configuration  of  Fig  5-1  will 
become  pseudo-singular,  whereas  the  second  configuration  constitutes  the  op¬ 
timum  choice  as  regards  the  orientation  of  the  conical  section  relative  to  the 

2 

computational  z-axis.  This  properly  is  demonstrated  by  the  added  sin  0  and 

2 

cos  0  -  curves  in  which  the  respective  =  sin  0  and  A  =  (cos  0  -cos  0  ) 

c  c  p  r  /u  r 

are  plotted. 

In  general,  the  optimum  distribution  of  the  aspect  angles  (0,0)  depends 

c  c 

upon  the  given  number  N  of  receiver  locations,  and  must  be  determined  for 

each  individual  N  separately.  Yet  from  the  general  structure  of  determinant 

C(0c>  N)  it  may  be  concluded  that  such  an  optimum  distribution  will  be  obtained 

if  the  aspect  angles  are  placed  in  the  now'  identical  maxima  of  the  Legendre  func- 
1 

tions  PN  (cos  0)  for  respective  N. 

5.4  COMPUTATIONAL  RESULTS 

With  the  evaluation  of  the  first  four  determinants  as  described  in 
section  3.2,  the  properties  of  matrix  C  (0  ,  N)  were  established,  where 

Det[c(0c,N))  =  Qn-  Pn(0c)  (5.4.1a) 

Qn  =  TT  Cc  (5.4.1b) 

c  =  1 
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Unit  Sphere 


Scatterer 


FIG. 5-l:"PSEUD0-SINGULAR"  DISTRIBUTION  OF  ASPECT  ANGLES 

(Invariant  axis  a  of  conical  section  (0  )  identical  with  the  computational 
z  -  axis ) 
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A 

Z 


COS  0 


FIG. 5-2: STABLE  DISTRIBUTION  OF  ASPECT  ANGLES. 

(Invariant  axis  a  of  conical  section  (6  )  oriented  perpendicularity  to  the 
backscatter  direction). 
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(5.4. lc) 


PM(0 

N  c 


N 

>=  TT 

c  =  1 


. 

sin  6 


N 

TT 

N>  r>s: 


(cos  0  -cos  6  ) 
r  s 


where  Cc  denotes  the  c  constant  multiplier  of  Eq.  (5.2.9a).  The  derivation 
of  a  general  algorithm  for  the  constant  multipliers  Cc  required  additional 
computation  of  the  next  higher  order  determinants  Det  C  ( 0  ,  N)  together  with 
the  corresponding  P^(9  ).  In  table  V-l  numerical  results  are  presented,  where 
Det  0(0^,  N)  was  computed  directly  from  the  matrix  formulation  as  derived 
from  (5.  1. 8a)  and  (5. 1 . 8b).  Apart  from  negligible  round-off  errors  the  computa¬ 
tional  results  verify  the  assumed  expression  of  C  to  its  best,  where 


_  [(2c  -l)'.'.]2 
C  [(c-1)'.]2 


2 

c 


N 


QN  ■  TT. 


c  =  1 


yy  c2  f(2c-  l)".]2 
c  =  1  [(c-1)'.]  2 


(5.4. 2a) 


(5.4. 2b) 


For  the  spherical  test  case  of  electric  measure  ka  =  2  (8579-3-Q)  the 

matrix  inversion  was  tested  for  a  varying  number  of  N  of  aspect  angles  and 

for  different  choices  of  the  sets  of  aspect  angles  (0  ,  0  ,  C  =  1,  2,  .  .  .,  N).  In 

c  c 

the  following  tables  V-2  and  V-3,  test  results  for  N  =  4  and  N  =  6  are  presented 
to  verify  properties  of  stability  as  they  can  be  predicted  from  the  determinant 
given  by  (5.  2.  9).  To  do  so  the  far  field  components  for  the  chosen  sets  of  aspect 
angles  are  computed  from  known  expansion  coefficients  as  given  by  Table  III-  of 
8579-3-Q  with  (5. 1. 6),  where  the  coefficient  extracted. 

Employing  the  matrix  inversion  technique,  as  described  above,  the  expan¬ 
sion  coefficients  are  computed.  For  a  computational  check  the  far-field  compo¬ 
nents  were  recomputed  with  these  coefficients.  It  was  observed  that  the  minimum 
deviation  occurs  for  those  chosen  sets  of  aspect  angles  (0c>  0c  =  45°,  C  =  1,  2, . . .  ,N) 
for  which  the  magnitude  of  the  determinant  is  of  the  order  of  unity. 
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COMPUTATIONAL  DERIVATION  OF  THE  CONSTANT  MULTIPLIERS 
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TABLE  V-2-la:  COMPUTATIONAL  RESULTS  FOR  N  =  4. 
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o 

0 . 5490327 

0.0 

-0. 7054T05 

0.0 

-1  .150044 

0.0 

4  0.0 

-0.24754080-05 

0.0 

-2.121316 

0.0 


57303291 


55 


TABLE  Yr2-lb;  COMPUTATIONAL  RESULTS  FOR  N=4. 


5  - 

0.4999998 

0.0 

0. R7665600-05 

o 

. 

o 

77575TJ7T 

0.0 

67TH7T95 

o.o 

6 

o.o. 

-0.7071080 

0.0 

-1.060658 

U.TT 

•  0  •  Z651  573 

0.0 

1 .1 04860 

7  -0.1  R30115 

0.0 

1.837124 

0.0 

ZTfHjHV) 

0.0 

-7.328176  — 

0.0 

8 

0.0 

-0.70fl0Rl 

0.0 

0. 74262520-05 

UTO  — 

1.060661 

u .  fJ 

-U . 1 856561 0-04 

~  A2  MATRIX  ..  ,  „ 

=  B(0  )  C 
c 

■  (0  )  (see  5.  2.  6) 

Tiro  .  . 

0. 6999979 

0.0 

■=0.87665280-05 

0.0 

-2.  62  5011 

0.0 

-6.187173 

Z 

u  rnrrrtT55 

“^TTD 

L  .060654 

0.0 

0.2651563 

0.0 

-1. 104R56 

0.0 

......  j 

“0.0 

0.1  8301  09 

0.0 

-1 . 8371 1  7 

o.o 

-Z.743R46 

0.0 

2.328167 

4 

“C. 7071 055 

0.0 

-  0  .  T4  Z6  ^  Z5TT-ir5  0.0 

-I . 060658 

0.0 

0.18565540-04 

0.  0 

0.0 

-o.  nrrro m 

070 

-1  .499998“ 

o.o 

-1 . 590980 

0.0 

-0.6249803 

~~E~‘ 

-0 .9555520'  “ 

“0.0 

1  .050558  ““ 

0.  0 

3. 844897 

0.0 

3.535501 

0.0 

7  ' 

“0.0 

-0.7071081 

8.  0 

“=0.5590347' 

0.^0 

0.7054130 

0.0 

1. 158048 

8 

0.  247541  m-06 

0.0 

- 2'.  121  324 - 

0.0 

-0. 40844340—04 

0.0 

-5.303310 

0.0 

DETA 

=  8.199224  = 

=  Det  {b  }  -  Det  [c  } 

(defined  by  5. 1. 

12) 

A  INVERSE 


1 

-24.18691 

43.94923 

-32. R3723 

10.47621 

“=74.95881 

54".  1  4295 

—7  7. 5  93  57 

9  .782475“ 

2 

24.95R91 

-44.14315 

?2. 59373 

-9.782481 

24  ,  1  8681 

-47.94904 

32. 83707 

-10.47616 

3 

P . 626^73 

-15.34219 

1  1  .  52027 

-3.555206 

7.  91  60  68 

-13.74816“ 

- 97759778 - 

“ =775T72T3“ 

4 

-7. 91 6098 

13. 74R22 

-9.  759324 

2.51 7228 

-R .  62659 4 

T5. 3421 2 

-11.52021 

T.  555185 

....  5 

-2.468427 

4.499455 

-3.373234 

1 • 006895 

-7. 463475 

4.137175 

-7.5V7613 

0.51 51431 

6 

2. 4634R5 

-4.137193  . 

2. 557626 

-0. 5151468 

Z. 46841 6 

-4.499433 

3.37371 1 

-1  . UO 68 89  “ 

7 

0.6154678 

-1.078557 

0.7737976 

-0.2207796  " 

0.4751979 

=075397716“ 

0 .7007681 

0715752550=05 

8 

-0. 4351990 

0.5392729 

-0.  2002686 

_ _  -  -  I _ -  .  _ !— _ 

-0. 16008140-05 

=075T5*5S5  7 1 .078551  =977737928““  0. 2257789 


TVF'T A 1  = - 8.557  744  =  Det  (B  }'»  'Det  fc^ 


Al FNVEPSE 

1 

8  •  T47T64 

0Y0 

=10.  54 7H7 

0.0 

-4 .1141  OR 

0.0 

-3.499777 

0.0 

T~ 

“7570“ 

0. 7642295 

070 

1.13 1351. 

0.0 

•  -3.016924 

0.0 

1 .979821 

=77372  744 

070 

11707329“ 

"0.0 

4.042033 

0.0 

3.254776 

0.0 

4 

“U.O  -1- 

-1.077481 

“0.0“ 

-1.279510 
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TABLE  V-2-lc:  COMPUTATIONAL  RESULTS  FOR  N  =  4 


0.0 

?. 693675 

0.0 

-3.232408 

b 

4  •  440(>3£ 

0.0 

-6. 500642 

~ nrrr” 

-2.298367 

0.0 

-1.687486 

0.0 

6 

”0.0 

0 .5028201 

~TT.1T 

- 0.7547380 

0.0 

-7.011284 

0.0 

2.262690 

7 

-T . 7  2  5368 

D  .IT 

T79TT5596 - 

0.0 

0. 8047444 

0.0 

0.4309940 

O.Q 

H 

"TT7TI 

-0.4309458 

070 

‘  -U.  32  324  71 

0%0 

1.077484 

0.0 

-1.292979 

DETA2  =  8.557728 

-  Det  B  •  Det 

°2 

A 2  INVERSE 

1  0.0  3.016933  0.0  -1.979826 

- 0.0 - =0.7542263 - 0.0 - -1 .131347 

2  4.U4134  0.0  3.499800  0.0 


-8.1471 57 

o.o 

10.54786 

0.0 

3 

0.0 

-7.693683 

0.0 

3.232417 

'  o  m 

L  •  07  f4  1 6 

rrrD 

1 .2  79506 

4 

-4.  042058 

0.0 

-3.754797 

0.0 

7.372738 

0.0 

-11.07378 

0 . 0 

5 

0.0 

2.011790 

0.0 

-2 . 262696 

~n*o  T“ 

-U. 507 El  TO 

0.0 

-0 . 7542349 

6 

2.  ?9'B32! 

0.0 

l »  68^497 

0.0 

-4 . 44U0  $  t 

u.u 

6. 50O436 

rr.n 

7 

c.o 

-1 .077487 

0.0 

1 . 292983 

o.o 

TT.  47  0OT4Q 

D .  0 

0 .3232456"  ~  ' 

8 

-0.8042488 

0.0 

-0. 4309969 

0.0 

1.275759“^ 

"0.U 

„  r  -l  1 

-1.905504 

0.0 

OE  I 

=  0 .95H 1 058 

*  DetlI  +  C!  C2j 

A  INVERSE 

l 

-?4. 1 840] 

/.a  .94923 

-32.83733 

10.47621 

-74.94881 

44.14295 

-32.503^7 

9.782475 

~T 

24.96491 

- 44 .14315 

377511371 

-O'.  782481  “ 

24.1 8681 

-42.94904 

32.83707 

-1 0.47616 

3 

8 • 42  65  73 

-15.34219 

'  11 . 52027 

-3 • 555206 

7.916048 

-13.74816 

9.759778 

-2.517213 

4 

=7191  6098 

13.74827 

-5.759824 

2.517228 

-8.626534 

15.34212 

-il. 52021 

3.555185 

5" 

—  2 . 4  4  M  4  7  f 

^“4.409455^ 

=373731 34 

1 . 00  67355 

-2.463475 

4. 1 37175 

-2.557613 

0.51 51431 

6 

~  7/463485 

^4“.TT7 1  93 

27357626 

“  -0. 5151468 

2.468416 

-4.499433 

3.373217 

-1.006889 

7  0.  6T546TR - -T.  078557 - 0.7737976 -  ^0.2207796  - 

0.4351979  -0.539271  6  0.7002681  0.  16752550-05 


'fl-0.'435  J.  990 - 0.5392  729 - -0.2002686 - -U.T600H1 40-06 

-0.6154646  1.078551  -0.7737928  0.2207780 


R  T  FR  ,  PI 

- - (Real  Part> - {imaginary  Parti 

1  1  0.47R?311  -0.3698460') 

7 - r  "0. 7822469 - -0.70409601 

9  1  1.075998  -0.1207969  ( 

"4  T  “I  .“772 995 - -071 5381 57T 

5  1  -0.2931354  0.5000026) 

6  I  — 074874  9T4 - 0732852737' 


E0 
—  M 


TABLE  V-2-ld:  COMPUTATIONAL  RESULTS  FOR  N  =  4  . 


7  I  -0.697716?  0.  1075299  )  Erf 

“R - 1  "-0.R6  7O1 47 - =0'."1  60Rq77'J —  VjJ 

-v  -  r  X9  •  xi 

(Real  Part)  (Imaginary  Part) 

“1  1  0.7.7?769R  -0. 00906]  ft 

2  1  0.6708078  0.4146107 

“7 - I~-n  .  7  099470 - ^Y5t77fiT74H=Tn 

4  1  ~0.  3776765  -0.2307197 

“5  T~7>.23B25?AF-or-0.9T4rt567F-03_l 

6  1  0.41  747 1 2  E~0 1  0  .300'’R57F-02  (  h 

"7 - 1  =0. 1  4  2  0  6  R 1  E—  OZ  ~  0.44973?0F-'05~T — ^elTT 

B  I  -0 .1  969776F-02  -0 . 86 1 7674F-05  J 


,  FRO 

--(Real  Part) 


1  0.4782111 


FIC 

-(Imaginary  Part) 
-0.14074*0 


2 

T~ 

0. 782/4  ft R 

-0.20409507 

rr1 

3 

1 

1 .07*097 

-0.1 207969  ) 

'  E0 

II 

4 

1 

L  .  V  l?W'~ 

-0 .157*158  J 

M 

5 

1 

-0. 29313*1 

0. 5000024  ^ 

ft 

I 

-rry;  4  8749  IT 

'  0.17R772 7~t 

,  Ew 

7 

1 

-0.6977160 

0.1035298  ( 

«u 
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TABLE  V-2-2:  COMPUTATIONAL  RESULTS  FOR  N  =  4. 


T 

R  THFTA 

PHI 

* 

1  " 

1  t. 0471  os 

0.785*979 

1 

2  2.094*94 

0.7853979  . 

l 

*  1.308995 

0.785*979 

1 

4  1.570794 

0.788*979 

OFTA  = 

410.7577 

=  Det  B  •  Det  C 

(see  5. 1. 12) 

DFTA  1  * 

1.075514 

=  Det  B  •  Det 

(see  5-.  2. 1  and  5.  2.  4) 

DFT/V2  = 

l .075574 

=  Det  B  «  Det  C? 

(see  5. 2. 1  and  5.  2.  6) 

OET  = 

38 1 .  $  1  36 

■ DetfI±^ 

•  (see  5.  2.  7b) 

T  F  R 

(ReaT'ParF)"' 

1  0.7R2?4*5 


PI 


(Imaginary  Part) 

-0.2040971? 


2 

3 

1  1  .0711  ft 1 

\  1  . 0759  97 

-0.594585R  T 

-0.1207976  f 

E 

e 

4 

5 

1  1.777995 

1  -0.4R74R74 

-0. 15**1 40  J 
0. *285*40  \ 

V 

6 

7 

1  -0.9*87797 

1  -0.(59*7119 

-0.7960107  r 

E 

0 

fU 

8  . 

1  -0.8679297 

-0.1695909  ) 

R 

(Real  V art) 

XI 

(Imaginary  Part) 

i 

? 

1  0.  7327468 

l  0.6708*06 

-0 ,°099419 
0.4146110 

' 

3 

4 

1  -0.2099591 

1  -0.3726835 

0.56766*25- 

-0.230**9* 

01 J 

ftoln 

5 

6 

1  0.238230*1-01 

1  0. 41  7^50VE-01 

-0.97^56645- 
0 .  *00  2  902  F— 

0*  s 
02 

u 

7 

8 

1  -0. 1420165E-02 
1  -0. 1970025F-02 

0.4491 2915- 
-0. R636775F- 

05 

•05 

b  . 

eln 

R 

T  E*C 

rir. 

1 

(Real  Part) 

1  0.78274*5 

(Imaginary  Part) 
-0.204Q975  \ 

2 

i 

7 

J.  1.071160 

1  T. 075904 

-0. 5945868  I 
-0.1207976 

L  Ee 

4 

5 

1  1.272994 

1  -0.487487* 

-0.1533140  „ 
0.3285239  > 

\ 

6 

7 

1  -0.9387795 

1  -0.693711 8 

-0.7960106  , 

0.1035333 

1 

X 

8 

1  -0.8679290 

-0.1695908  v 

) 

M  - 
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TABLE  V-2-3:  COMPUTATIONAL  RESULTS  FOR  N  =  4. 


T 

P 

TMFTA 

PHI 

1 

} 

1.178097 

0.7R53979 

1 

2 

1  .308905 

0.7853979 

1 

3 

1  .439R96 

0.7853979 

] 

4 

l  .  570794 

0.7853979 

OF.T'I  = 

0.  1 172266F-0 1 

=  Det  B  <  Det  C 

(see  5. 1.12) 

=  Det  B  •  Det  C 

(see  5. 2. 1  and  5. 

DFTA1 

=  0. 

1 375974 

OFTA? 

=  0. 

1 3759RR 

=  Det  B<  Det  C0 

(see  5. ‘2. 1  and  5. 

OFT  = 

0 . 8  5 1 9  5 1 8F-01 

-  De‘{i  +  Cil '  C21 

(see  5. 2.  7b) 

R 

T 

T9 

5  I 

(Real  Part) 

(Imaginary  Pa3*t) 

i 

i 

0.036435? 

-0. 1 500959  1 

2 

i 

] .075997 

-0.1 207976  /  F 

3 

i 

] . ) 9241 0 

-0.1 210306  (  9 

4 

i 

1.272995 

-0.1538140  ) 

5 

i 

-0.591 56^0 

0.2225537  \ 

6 

i 

-0.69?71 19 

0.1  O'1  5  RR  3  ( 

s  E  w 

7 

i 

—0 . 7R7  97  01 

— 0 .2 76 5677 F— 01 

% _ 

8 

i 

-0. 8679292 

-0 • 1 6OS909  ^ 

/  1 

XR  '  XI 

(Real  Part) _ (Imaginary  Part) 

0.7**?  74  0  8  -O.RC°o416 

0.6708390 _ Q.4  1  4MJ  5 

- 0  .  20Q95  5 9  0. 5676645»=-0  1 

-0.  3776856 _  -Q.  ??07R9l _ 

0.2387?3*f:-01  -0 *9745671 T- 03  3 
0.4  1749969-01  0.300? 943F-02  (  b 

-0. 1 4i98?0F-02  0  •  45  6092 1 F- 05  (  eln 

-0.19701  8 1 F  -  0?  -0.8590585F-Q5  , 


R _ _ T _  FRC 

.  (Real  Part) 
1  1  0.9354351 

?  .  1  1.075996 

3  1  1.1 9241 0 

4  1  1.272994 

5  1  -0.5915629 

6  1  -0.6937119 

7  _ 1  -0.7879701 

8  1  -0. R679291 


If. 


(Imaginary  Part) 
-0.1 90095  Q 


-0.1 207Q76 
-0 .1 21 9306 


-0. 1539140 

0.2225536 


-I u- 


0.1035333 
-0. 2  7656  77F-01 


LfcL 


-0.1695908 
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TABLE  V-2-4:  COMPUTATIONAL  RESULTS  FOR  N  =  4. 


T 

R  THFTA 

PHI 

1  ' 

1  0.5235986  . 

0.7853979 

1 

2  0.7863979 

0.7853979 

i 

3  1.047196 

0. 7853979 

i 

4  1.308995 

0.7853979 

oeta  = 

0.3053473 

0ETA1  = 

0.  M  43691 

DETA2  = 

0.6143843 

OET  =  0 

.4970094 

P 

T  FR 

FI 

1  . 

1  0.2273?  '*4 

-0.631  *-700 

2 

1  0.47B?7°R 

-0.3603482 

3 

1  0.78224^5 

-0. 20'.0O76 

u 

1  1.0769O? 

-0.1 207976 

5  . 

1  -O.n^oiP1 

0.61 P605R 

6 

1  -0. 2931 336 

0.5000023 

7 

1  -0.4874874 

0.3286240 

8 

1  -0*6977119 

0.1036333 

R 

T  XP 

XI 

i 

1  0.7327479 

-0  .Q099606 

2 

1  0.6708288 

0.41^4? 03 

3 

l.  — 0 . 2  0°95  o  i 

0.567  H  89E-01 

4 

1  — 0 .  3"/26  82  5 

-0.2303446 

5 

1  0. 2382298F-01 

-0.97S9669F-03 

6 

1  0.417445RF-01 

0.30041.37E-02 

7 

1  -0. 14701 20F-02 

0 .4689246^-05 

R 

1  -0.196°866F-02 

-0. 87P^497F-05 

R 

T  FR  r. 

FIC 

1 

1  0.2273245 

-0.6716700 

2 

1  0.47 82 296 

-0.3 60 74 8 2 

3 

1  0.7822435 

-0. 704Q9 76 

4 

1  1 ,075906 

-0.1207977 

5 

1  -0.17991R7 

0.6186067 

6- 

1  -0.2931336 

0.5000023 

7 

1  -0.4874873 

0.3285239 

8 

"l  -0.69371  18 

0.1035333 
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TABLE  V-3-la:  COMPUTATIONAL  RESULTS  FOR  N  =  6. 


t — ~  '“R - rHFr/r- - imT 


1 

l 

T 

1 

T 

1 


"1 0.78540  00* - 0.7894000 

2  1.047!°°  0.7854000  NR 

*3 - 17708999  0.7854000 

4  1.570800  0.7854000 

“5"  '  0.? 0770750 - 0  .7 0770 (TO - 

4  0. 5214000  0.7854000 


DF  *  0  NR  =  6  NT  *  1  Computation  of  Matric  Elements 


T  -  18=  1  THFT  A ( T  *  R )  =  0.7854000  PH1(T,RI  =  0.7854000 

- 1 - 0.7071055 - 07707178 1 - 070 - 

2  0.7409957  1.499999  1.500005 

- 1— -0.  1757801  1.590981  5.703708 - - 

4  -0.406251 2  0.4249819  9.774077 

- 5"“-0.  7756481  -  *.-0.0943888  I  ~“  9.380475 - -  - — 

4  -0.1484324  -2.296881  1.440434 

7  5 - T7F*= - 7  1 H  FT  ATT  7R 1  = —  1.04  7199 - PHTTTTR1  = - 077854000 - 

1  0.4999985  0.8640243  0.0 

- 2— -0.1250027  ~*  1.299035  0.750004 - .  : 

3  -0.4375007  0.3247495  5.62409] 

- S- -0.  2890506  -  -1.7531  77  “  4.21  8695  ~*  <  - 

5  0. 89847*3F-01  -1.928251  -4.921961 

- 4 - 0.723  2  42  8 - =77  4 07 751 *8 - =74. 15070 - 

T  =  1  R  =  7  THFT  A(  T  .8 1  =  l. 308009  »HHT,R|  =  0.7854000 

- 1 — 0.2588170  0.9650965  0.0  . .  ~  - 

2  -0.3995207  0.749Q944  2.790941 

- 7~0. 7448823  - -0.0636111  3.422195 

4  0.1424370  -1.581923  -1.71  4404 

5  077577750  077777307  -10. 12.999 

6  0. 4309556F-01  2.059417  -0.7062795 

7  5  1  RV  4  THE7A(T,R)  =  1.570800  PHT(T,R7“a - 077854000 - 

1  -0. 3500762F-05  1.000000  0.0 

- ?~0. 5000000  “*  ~ =0.105022  98-04  “  3.000000  - 

3  0.52511  40F-05  -1  .500000  -0 . 5 2 5 l 1 40F- 04 

- 4 - 073750000 - 0.26755708=04 - =77700700 - 

5  -0.65639238-05  1.875000  0. 1 837899E- 03 

- 6~0.*3125000 - -0 . 4  5  947488-05 -  13.12500 - 

T- »  1  R  a  5  THFT A( T 1 R 1  =  0.7927000  PHI (T  ,R1  »  0.7854000 

- 1 - 0.9238703 - ; — 0.3826847 - 0.0  - 

2  .0.7807288  1.060661  0.4393414 

- 3 U.  68562  9  8 - 17137 5781 - “7..  02949] - : - 

4  0.361 597  R  2.629439  5.464164 

5  0.T3751  86 - 3.1  21 167*  "-= - 17  .0356° - 

6  -0. 7636195F-01  3.188541  18.60277 

1  * - r  R*  a  F  THFTATT7R  1  =  0.5736000  PHI TT »R1  a - 077854000"" 

1  0.8660249  0.5000010  0.0 

2  0767*59781  l .  799TTTg  077500079  ~ 

3  0.3247566  2.062498  3.247604 

4  - 0.234T4T6E-01 — “2.435680 - 77* 96R75H - 

5  -0.2237749  2.167948  14.20321 
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8 

1 

-0. 1962313E-02 

-0. 1 2877R7F-04  ^ 

9 

L 

0T5TRBT  H6fc-u4 

-U  •  L  0 1 4 1 00 t— Ob  f 

10 

1 

0.63R6307E-04 

0.92  7985  8F-06  } 

b  n 
el  n 

1 1 

1 

~0'.T5*n  RFRF-7)5“ 

0.12781478-06  ( 

12 

1 

-0.141 1 963F-05 

— 0.1046130F— 06  / 

R 

T 

FR  C 

FIC 

1 

1 

0.' 

2 

1 

0  . 

3 

1 

1 

4 

1 

1 

5 

1 

0. 

6 

~r 

u. 

7 

i 

-0. 

R 

V 

~-0  . 

9 

l 

-0. 

1  0 

'  i 

-0. 

11 

l 

-0. 

■0.3(S0TA(S0 
•0.  7040949“ 
-0.1  207967 

■0.133  81  57“ 

*0 . 604  fl  7  7  7 
=DT5TT7T6W 

0.3000025 
~0;i7R57?3~ 
0 . 10*5299 
-0.1  693  077' 
0.6589078 
"0'.6TH50'64- 


TT 


T ' -- 0 .13  974 9T 


68 


TABLE  V-3 

-2: 

COMPUTATIONAL  RESULTS  FOR 

T 

P 

THFTA 

PHI 

* 

1 

1 

1 . 570794  • 

0.7853979 

i 

? 

1 .400096 

0.78R3979 

i 

3 

1  .  30R9n* 

0  •  7  R  *  '*  9  7  9 

i 

4 

1  .1 7809  l 

0. 7R53979 

1 

5 

J .0471 06 

0.7863979 

i 

6 

0.°162976 

0.7863979 

Dp  TA  ^ 

-0.24 

Ri  404F-0’ 

=  Det  B  Det  C 

=  Det  B  Det  C 

Df^TAl 

=  -0  *  946 1 M IF- 07 

1 

HETA2 

=  -  0  *  9 

46!  81  6F  -  0  * 

=  Det  B  Det  C( 

2 

OFT  = 

0.454335* 

=  Det  |l  +  C 1  C2j 

P 

T 

-  FR 

C  I 

(Real  Part) 

(Imaginary  Part) 

1 

1 

1  .777Q70 

-0  •  1  c  7  8 1  4  0  ' 

V 

2 

1 

1 . 1 027  7  7 

-0 .121030s 

3 

i 

1  . 0764  77 

-0.1 707Q75 

V  Fj  -j. 

4 

1 

0.92*8061 

-0.1 5009*9 

r® 

6 

1 

0. 78737?o 

-0. 20 4097 6 

6 

1 

0 .67667  /, fl 

-0.2765647  , 

) 

7 

1 

-0. R6R0?  RO 

-0.1 69591 0 

8 

] 

-0 .7RR?RK0 

-0.276567*=- 

-01  / .. 

9 

1 

-0. 6040970 

0.1076372 

 }-** 

1  0 

i 

-0. *91 8370 

0.2226577 

1  1 

i 

-0.48761?? 

0.3285241 

 1  . . .   

12 

\ 

-0. 3863768 

0.4210207 

) 

P 

T 

XR 

XI 

(Real  Part) 

(Imaginary  Part) 

1 

1 

0.7377083 

-0. Q1 007S? 

\ 

2 

1 

0.6708680 

0.4146963 

j 

3 

1 

-0 . 20004?/. 

0.5679631= 

— 0 1  (  a 

4 

1 

-0.377693R 

-0.  ?  507743 

>  oln 

5 

_ 1 

0.23816  74F- 

01 

-0.98647436 

-0  3  T 

6 

1 

0.41  747  P6p  — 

01 

0.301 49676 

“0?  ) 

7 

1 

-0. 141751 9F- 

02 

0.7196473= 

-05  S 

8 

1 

-0.1 t'6°6R6F- 

03 

-0.1422371  f= 

-04  / 

Q 

1 

0. 5137791F- 

04 

-0.6797777= 

-06  >  b  , 

1  0 

1 

0. 65 iO°74E- 

04 

0.  1 60 9 R ]  Rr 

-05  (  eln 

n  1  -0.1117?OOF-O5 

12  1-0.  13821  73E-05 


■0.581  '’6086-07 
-0,37211 85F-0 6 


R 

T 

fr  r 

8  1  C 

(Real  Part ) 

(imaginary  Part) 

1 

1 

1  .2730  7R 

-0.15*3140  ) 

2 

1 

l  .192776 

-0.121  '70  6  T~ 

3 

1 

1 .076476 

-0.1207975  V 

0 

4 

1 

0.9*58060 

-0.1500969  r 

oln 

5 

1 

0.782322° 

— 0. 7049976  1 

6 

1 

0.6266248 

-0.7765667  ) 

7  . 

1 

-0.8680279 

-0.169=910  N 

8 

1 

-0.78H284R 

-0.27  65 6765-01  I 

b  , 

9 

1 

-0.6940969 

0.1035373  / 

eln 

1  0 

1 

-0.591  8369' 

0.2225537 

69 


TABLE  V-3-3:  COMPUTATIONAL  RESULTS  FOR  N  =  6 


T 

p 

thfta 

PHI 

1 

1 

o.i ?nagos 

0.7847979 

i 

? 

0.?61 7993 

0.7868979 

1 

3 

0.  3026989 

0.78*3979 

1 

4 

0.5?3*9P6 

0.7853979 

■ - 1 — 

6 

0.654498? 

0.7853979 

i 

6 

0.78*3970, 

0.7863979 

HFTA  = 

-0.1  06076  9F-1  (S 

DFTA-1  * 

-0. ??40667F- 1 5 

OFT  A?  = 

-0 .6206R44F-1  6 

1ZT  = 

0. *77847 3F-01 

P 

T 

PR 

FI 

l 

1 

1 

0.7  564  7  95  F- 01 

-0. 607^904 

? 

1 

n.666?36.F-01 

-0.( 61 6657 

3 

1 

0.1 347047 

—  0 . ( 0487Q6 

4 

i 

0.??701 Ri 

-0.531670? 

5 

1 

0. 3477306 

-0.4476709 

6 

1 

0.47771 06 

-0.36034^3 

7 

\ 

-0 .2000466E-01 

0.70^730 

8  « 

) 

-0.441  07]  1  *=-0l 

0.6*75667 

9 

i 

-0.B427298F-01 

0. 6 6  8 r< 0^6 

1  0 

i 

-0 • 1307481  ’ 

0.6185058 

11 

i 

-0.2097238 

0. 565^968 

1  ? 

1 

-0. 797698° 

0 • 600  00?  3 

p 

T 

XR 

XI 

1 

1 

0. 753*4? 8 

-0. 700*?06 

? 

1 

0. 65?RR?4 ' 

0. 7Q76AQ3 

3 

1 

-0.71  77883 

-0.1 08*31 3r-0 1 

4 

1 

-0.36B4531 

-0.1 630370 

_  5 

1 

0.26601 P0F-01 

0. 1 9P85R7F-0! 

6  , 

1 

0.41  068*85-, )] 

-0 . 1  8  76  36 3*- 01 

7 

1 

-0.71 799R?F-02 

-0.61?4780l=-0? 

R 

l 

-0 • 70606  74  F  — 0  3 

0. 5146421 g-n? 

Q 

i 

0. 19H°400r-03 

0.81078089“ 03 

1  0 

i 

0. 1 207975F-0? 

-0.8477415F-03 

11 

1 

-0. 1 48852OF-04 

-0. 65°02?9It-04 

1? 

1 

-0.9O494O3F-05 

0.68161 87c- 04 

R 

T 

FRC 

nr 

1 

1 

0. 748O265E-01 

-0.66R?! 1? 

2 

1 

0 .66RC2 685-01 

-0.661 9473 

1 

0. 1 3361 78 

-0.4046321 

4 

1 

0 . 7766036 

-0.5307534 

5 

1 

0. 347401 1 

-0  .4A6??1  1 

6 

1 

0.4775949 

-0.3579487 

7 

1 

-0.191 19O7F-01 

0.70^*9143 

a 

1 

43?08RRF-01 

0 .6P77364 

Q 

i 

-0. 833P946E-01 

0.6589058 

1  0 

i 

-0.1388899 

0*6181 347 

70 


TABLE  V-3-4:  COMPUTATIONAL  RESULTS  FOR  N  =  6; 


T 

P 

THETA 

PHI 

1 

i 

1  .  67  07  9  4 

0.788? 979 

i 

2 

1  .308995 

0.7853979 

3 

l. 832505 

0.7853979 

i 

4 

1.0471  Of; 

0.7883979 

1 

*5 

?. 094794 

0 .7B51979 

1 

6 

0 • 7  85?  9  70 

0.785^979 

OT.TA  = 

-l  766898. 

oetai  = 

-0. 1 6479flflr-07 

DFTA2  = 

-0.164A334E-02 

OFT  * 

0.10721  8  3  F  10 

R 

T 

TP 

C1 

1 

1 

1  .273070 

-9.  1  *3R1  4  0 

2 

1  • 

1  •  07f>4  77 

-0 • 1 307976 

i 

1 

1  .286344 

-0. 31 P 3091 

4 

1 

0 . 7823720  • 

-0.2040975 

5 

l 

•  1  .  0  7  0888 

-0. 594  5  84  8 

6 

1 

0.4 777 1  OS 

-0. 06O’4R3 

7 

1 

-0.R6B0-2B0 

-0.1  69ioj  0 

P 

1 

-0. 6940970 

0. 10^4333 

9 

1 

-o. 96nssss 

-0.4770909 

1  0 

1 

*0 .48751?! 

0.3?Rr<741 

1  l 

1 

-0. 03852  33 

-0.796^1 07 

1  2 

1 

-0,  ?0 ?/,98? 

0.4000<r,1 

R 

T 

XR 

XI 

1 

1 

0 . 777741 7 

-0.90004 1 6 

7 

i 

0. 670R1  7S 

0  .  M  '■  4 1  0  9 

3 

i 

-0.  ? 0 *3 0 5 6 0 

0.567mA  •»7P-01 

4 

1 

-0.3 ^8 58 

-0 .7307191 

5 

1 

o.73R2i PSF-cr 

-0 .974  '.91  1  P-01 

6. 

1 

0.4174646F-01 

0. 3007P73E-0? 

7 

i 

-0. 141 94185-0? 

0.4501 417P-05 

8 

1 

-0. 19703 74 r-0? 

-0 . RS  n07 2 9 r-0S 

9 

l 

0.S178S7-E-04 

0. 45805605-08 

1  0 

l 

0.656861 7F-04 

0.1160489^-07 

1 1 

1 

-0.  1 2076?  7F-05 

0.666R1 04F-0B 

12 

1 

-0.1 S6R01 7E-0S 

0 . 1  0?  7  6r>  flE-07 

R 

T 

FRC 

PTC 

i 

1 

1.  .  ?  7?  0  7R 

-0. 1538140 

7 

1 

•  3.  .0  764  74 

-0.  1.207976 

7 

•i 

•  1.286345 

-0.3183992 

4 

1 

0. 7873?  ?o 

-0.2040974 

5 

1 

1 .070887 

-0 .  5  9^  5  R4  R 

* 

1 

0.4  7771  OS 

-0.360  A4R3 

7 

1 

-0.  P6S0?  70 

-0.1695909 

8 

1 

-0.6940970 

0.1025337 

9 

'  1 

-0.9605554 

-0. 477090B 

10 

1 

-0.48751 21 

0.32B5240 

71 


TABLE  V-3-5:  COMPUTATIONAL  RESULTS  FOR  N  =  6 


T 

P 

TMFT  A 

PHI 

i 

i 

l  .1  78  0 9  7 

0.7860979 

i 

7 

'  ,04  71  os 

0.7040979 

1 

7 

0.0167974 

0 . 7850979 

i 

4 

0. 7883079 

0.78*7979 

J- 

S 

0.6 4440  BO 

0.7860979 

i 

6 

0.6235986 

0.7850979 

ofta  * 

-0. 1 0  3  0  H  ?  0  F  -  0  * 

0FT41.  = 

-0.9475437F-06 

OFTA?  = 

-0. 9465356F-06 

OFT  = 

0 . 1  08  78  9^ 

P 

T 

FR 

«-  r 

1 

1 

3.0348061 

-0. 1  ^ 00° 8  0 

? 

1 

0.78^3779 

-n .  7940975 

3 

0. 6?  4  6? '•  8 

-0.0764667 

A 

1 

0.4777108 

-0. 0600480 

5 

1 

0.  V-273P5 

-0.4479700 

6 

1 

9. 2270  1  pi. 

-0.531 *70? 

7 

1 

-0. 891  8770 

0. ???^77 

8 

1 

-0.4R741 ?7 

0 . 30P  6041 

9 

1 

-0. 3  P  40  7  4  M 

0.4219207 

1  0 

1 

-0. 

0. 4090003 

1  1 

1 

-0 • 709  77  7  8 

0. 5654068 

\? 

1 

-0. 1397482 

0.61  88Q8  3 

R 

T 

XD 

XI 

1 

1 

0.70] 4806 

-0. 9 1 ??691 

? 

1 

0.470041 ’ 

0.41 7008Q 

7 

1 

-0.7098778 

0.5746462P-91 

4 

1 

-0. 3731167 

-0. oo i i 934 

5 

1 

0.23700 79F-01 

-0.  j.  21  Q506E-00 

6 

1 

0.41 B06 04 P-01 

0. 0300967E-0? 

7 

1 

-9. ] 0921 3  4  F - 0  2 

0. 613790  6  OF- 04 

8 

1 

-0 . 20i O424E-02 

-0. 93  5  6  30  8  F- 04 

9 

1 

0. 47R04PRE-04 

-0.926981 7F-Q5 

10 

1 

0. 741 61 7\ P-04 

0.17186R5E-04 

i  ] 

1 

-0.10  704  60F-05 

0.61 0^7768-06 

1  2 

1 

-0. 25034C1 F- 05 

-0. 1 88  M  1  ir-05 

R 

T 

FRC 

ric 

i 

1 

0.  <->388058 

-0. l 500087 

? 

1 

0.  78^71', 

-0.2040980 

3 

1 

n. 6044046 

— 0 • 7764674 

4 

l 

0.4  7771  14 

-0. T 60^4 9? 

S 

1 

0.7497404 

-0.447871 7 

6 

1 

o. 2270 l on 

-0.57  1  A  7  09 

7 

1 

-0.5918363 

0. ??7r 530 

8 

i 

-0.44741 16 

0.3?85?73 

9 

i 

-0. 0P607 64 

0. 4210700 

10* 

1 

-0.P9769R1 

0.500001 8 

72 


REFERENCES 


Barrar,  R.B.  and  C.L.  Dolph  (1954)  "On  a  Three  Dimensional  Problem  of 
Electromagnetic  Theory,"  Journ.  Rat,.  Mech.  and  Anal.,  Vol.  3,  No.  6, 
725-743. 

Leontovich,  M.A.  (1948)  Investigations  of  Propagation  of  Radio  Waves,  Part  II, 
Moscow. 

Morse,  P.M.  and  H.  Feshbach  (1953)  Methods  of  Theoretical  Physics,  Vol.  II, 
McGraw-Hill  Book  Co. ,  New  York. 

Muller,  C.  (1956)  "Electromagnetic  Radiation  Pattern  and  Sources,  " 
IRE-TR-AP-4,  No.  3,  224-232. 

Senior,  T.  B.A.  (1961)  "Impedance  Boundary  Conditions  for  Imperfectly  Conduct¬ 
ing  Surfaces,  "  A££b_Scii_ResiJ__Section_B,  Vol.  8,  418-436. 

Senior,  T.B.A.  (1962)  "A  Note  on  Impedance  Boundary  Conditions,"  Am. 

Journ.  of  Physics,  Vol  40,  663-665. 

Stratton,  J.A.  (1941)  Electromagnetic  Theory,  McGraw-Hill  Co. ,  New  York. 

Weston,  V.H.  (1963)  "Theory  of  Absorbers  in  Scattering  Theory,  " 

IEEE-TR-AP,  Vol.  11,  No,  5,  578-584. 

Weston,  V.H.,  J.J.  Bowman  and  Ergun  Ar  (1966)  "Inverse  Scattering  Investi¬ 
gation  -  Final  Report,  1  Oct.  1965  -  30  Sept.  1966,  "  The  University  of 
Michigan  Radiation  Laboratory  Report  No.  7644-1-F,  AF  l9(628)-4884. 

Weston,  V.H.  (1967)  "Inverse  Scattering  Investigation  -  Quarterly  Report 

No.  1,  April  1966,  "  The  University  of  Michigan  Radiation  Laboratory  Report 
No.  8579-1-Q,  ESD-TR-67-517,  Vol.  I. 

Weston,  V.H.  (1967)  "Inverse  Scattering  Investigation  -  Quarterly  Report 

No.  2,  July  1967,  "  The  University  of  Michigan  Radiation  Laboratory  Report 
No.  8579-2-Q,  ESD-TR-67-517,  Vol.  II. 

Weston,  V.H.,  W.M.  Boerner  and  C.  L.  Dolph  ( 1967)  "Inverse  Scattering 
Investigation  -  Quarterly  Report  No.  3,  Nov.  1967,"  The  University  of 
Michigan  Radiation  Laboratory  Report  No.  8579-3-Q,  ESD-TR-67-517, 

Vol.  IH. 

Weston,  V.H.,  W.M.  Boerner  and  D.R.  Hodgins  (1968)  "Inverse  Scattering 
Investigation  -  Quarterly  Report  No.  4,  Feb.  1968,  "  The  University  of 
Michigan  Radiation  Laboratory  Report  No.  8579-4-Q,  ESD-TR-67-517, 

Vol.  IV. 

Weston,  V.H.  "Theory  of  Absorbers  in  Scattering’,'  Report  CAA-0020-10-TR, 
Conductron  Corporation,  Ann  Arbor,  Michigan. 

Wilcox,  C.H.  (1956)  "An  Expansion  Theorem  for  Electromagnetic  Fields," 
Comm.  Pure  Appl.  Math.,  9,  115-134. 


73 


DISTRIBUTION 


Electronic  Systems  Division 

Attn:  ESSXS  27  copies 

L.G.  Hanscom  Field 
Bedford,  Mass.  01730 


Electronic  Systems  Division 

ESTI  23  copies 

L.G.  Hanscom  Field 
Bedford,  Mass.  01730 


Unclassified 

Security  Classification 


DOCUMENT  CONTROL  DATA  -  R  &  D 


1.  ORIGINATING  ACTIVITY  (Corporate  author) 

The  University  of  Michigan  Radiation  Laboratory,  Dept,  of 
Electrical  Engineering,  201  Catherine  Street, 

Ann  Arbor,  Michigan  48108 

2a.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

26.  GROUP 

N/A 

3.  REPORT  TITLE 

INVERSE  SCATTERING  INVESTIGATION 

FINAL  REPORT 

4.  DESCRIPTIVE  NOTES  (Type  ol  report  and  inclusive  dates) 

Final  Report  (3  March  1967  to  3  March  J968) 

5  Au  THORISff^irsf  nama,  middle  initial,  last  name) 

Vaughan  H.  Weston 

Wolfgang  M.  Boerner 

6.  REPOR  T  D  A  TE 

April  1968 

7a.  TOTAL  NO.  OF  PAGES  7b.  NO  OF  REFS 

73  14 

8a.  CONTRACT  OR  GRANT  NO. 

FI9628-67-C-0I90 

b.  P  ROJ  EC  T  NO. 

C. 

d. 

8a  ORIGINATOR'S  REPORT  NUMBER(S) 

ESD-TR-68-215 

8b.  OTHER  REPORT  n  0(5)  (Any  other  numbers  that  may  be  as  sign  ad 
this  report) 

8579-1-F 

10-  DISTRIBUTION  STATEMENT 

This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  unlimited. 

11.  SUPPLEMENTARY  NOTES 

12.  SPONSORING  MILITARY  ACTIVITY 

Deputy  for  Surveillance  and  Control  Systems, 
Electronic  Systems  Division,  AFSC,  USAF, 

LG  Hanscom  Field,  Bedford,  Mass.  01730 

13.  ABSTRACT 


A  brief  review  of  the  salient  features  of  the  theoretical  investigation  of  the  c.  w. 
bistatic  inverse  scattering  problem  is  presented.  The  effect  of  changing  the  origin 
of  the  coordinate  system  upon  the  convergent  properties  etc. ,  of  the  spherical  vector 
wave  function  representation  of  the  near  scattered  field  and  the  surface  loci  Ex  E_  =0, 

is  discussed.  It  is  pointed  out  that  a  great  deal  of  analysis  remains  to  be  done  in  this 
area.  The  determination  of  the  surface  of  the  scattering  body  from  knowledge  of 
the  local  total  electric  field  is  given.  Emphasis  is  placed  upon  the  generalization  of 
the  condition  Ex  E  ‘=0  as  applied  to  perfectly-conducting  bodies,  to  scattering 
surfaces  characterized  by  the  impedance  boundary  condition.  Properties  of  the 
matrix  inversion  associated  with  the  determination  of  the  expansion  coefficients 
from  far  field  data  are  discussed.  Some  numerical  results  are  presented,  and 
restrictions  upon  the  choice  of  aspect  angles  are  deduced. 
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